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Abstract 


In this work the formulation of structural dynamic equation of motion for a 
general helicopter rotor blade with swept tips is attempted. The work aims to solve the 
problem of the bearingless and the hingeless configuration of the rotor blade. The 
equations of motion have been derived using Hamilton’s principle and the Finite Element 
Method is applied for its solution. The formulation is validated by comparing the results 
of the present analysis for a uniform hingeless rotor blade with that of those available in 
literature. Results have also been generated for practical helicopter rotor blade. 

The effect of root offset on the natural frequency and mode shapes of the 
hingeless rotor blade has also been analysed. The coupling effect of geometric pitch on 
the natural frequency and mode shapes of the practical rotor blade has also been 
analysed. 
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Chapter 1 
Introduction 


The time varying loads on the main rotor system contributes significantly to the 
vibration in helicopter. Therefore, the structural dynamic characteristics of the rotor blade 
and also the dynamic characteristics of the fuselage have a very strong influence on the 
vibratory levels in helicopter. Any analytical study pertaining to the dynamics of 
helicopters requires the development of suitable mathematical models for: 

• Rotor blade 

• Fuselage 

• Rotor-fuselage interface 

Rotor blade model essentially consists of the development of structural, inertial 
and aerodynamic operators associated with its motion. The fuselage model is represented 
by an idealized structural model of a three dimensional structure. The rotor-fuselage 
interface model must represent both the geometry of the interface as well as the 
aerodynamic interaction in an appropriate manner. The fundamental tool, necessary for 
the analysis and design of helicopters, is the rotor blade dynamic model. 

During operation, the blades experience large bending and centrifugal loads. In 
order to relieve the rotor bending moments experienced by the blades, early rotor blades 
were provided with flap and lag hinges at the root of the blade. In addition, a pitch 
control bearing was provided to control rotors. Such rotor systems are usually referred to 
as articulated rotors. A schematic diagram of an articulated rotor system is shown in 
Fig. 1.1. The large number of moving parts leads to a mechanically complex rotor hub 
system accompanied by the associated wear out problem requiring fi-equent maintenance 
and replacement of parts. With advancement in technology, increasing emphasis has been 
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placed on the development of hingeless rotor systems The construction of these rotors 
relatively simple, because of the absence of flap and lag hinges; but a pitch bearing is still 
provided for blade pitch control A schematic diagram of hingeless rotor system is shown 
in Fig 1.2 The hub assembly and the main rotor blade geometry of practical helicopter 
are shown in Fig 1 3. 

1.1 Structural Modeling 

Since the rotor blades are long, slender beams, they will undergo moderate 
deformation. A non-linear strain-displacement model is used to describe the coupling 
effects between axial, bending and torsional deformations. Generally, the strains are 
assumed to be small in comparison to unity Such an assumption is consistent with the 
design requirement based on fatigue life consideration which states the rotor blades must 
be designed to have an operating strain level well below the elastic limit of the blade 
material 

In order to develop an indigenous rotor blade dynamic model, Punit Kumar Gupta 
[1] formulated a finite element of a hingeless rotor blade under-going flap, lag, torsion 
and axial deformations. The rotor blade was treated as a straight uniform beam. 
Subsequently, Venu Gopal [2] extended the model by including non-uniform properties 
of the blade and swept tip These rotor models included all the complex geometric 
parameters such as torque offset, blade root offset, pre-cone angle, pre-droop angle, pre- 
sweep angle, pre-twist, tip sweep angle, tip anhedral angle and hub motion. However, in 
these studies rotor speed was taken as constant. 

During starting up, the helicopter rotor experiences variable speeds. In order to 
consider transient conditions in rotor r p.m, it has been treated as a variable. 
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1.2 Objectives 


Objectives of present study are 

1 To develop a most general beam type finite element model for helicopter 
rotor blade, including all the complex geometric parameters such as torque 
offset, blade root offset, pre-cone angle, pre-droop angle, pre-sweep angle, 
pre-twist, tip sweep angle, tip anhedral angle and hub motion 

2 To validate the model by comparing the results of this study with that of 
those available in literature. 

3. To check the model for practical data, provided by industry. 

4. To conduct detailed studies on the dynamic behavior of rotor blade to 
determine the effects of root offset and geometric pitch on natural 
frequencies and mode shapes 
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Chapter 2 

Rotor Blade Model and Assumptions 


Helicopter rotor blades are long slender beams attached to the hub through a complex 
geometrical and mechanical arrangement. The geometrical parameters describing the 
configuration of the rotor blade-hub system is shown in Fig. 2.1. The parameter a 
represents the torque offset, which is the distance from the center of rotation (hub center) 
to reference axis of the blade. ei and ez represent blade root offset distances from the hub 
Pp represents the pre-cone angle defining the orientation of the blade pitch control axis 
with respect to the hub plane. P* and Pd represent the pre-sweep and pre-droop angles, 
respectively, and representing the inclination of the blade reference axis with respect to 
the pitch control axis. The blade consists of a straight portion and a swept tip whose 
orientation relative to the straight portion is described by a sweep angle As and an 
anhedral angle Aa 

2.1 Basic Assumptions 

In the formulation of the dynamic model of the rotor blade with swept tip, following 
assumptions are made: 

1 . The blade is treated as an elastic beam. 

2. The blade is modeled by beam finite elements along the elastic axis of the 
blade 

3 . The rotor shaft is rigid 

4. The blade cross section has a general shape with distinct shear center, 
aerodynamic center and center of mass. 

5 The blade undergoes moderate deformation in flap, lag, torsion and axial 
modes. 

6. The blade has non-uniform properties along the span though it is made of 
isotropic material. 
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2.2 Ordering Scheme 


In the formulation of the equations of motion of a rotor blade with swept tip 
undergoing moderate deformations, a large number of higher order terms are generated 
In order to identify and eliminate higher order terms in a consistent manner, an ordering 
scheme is employed. This ordering scheme is based on the assumption that the slopes of 
the deformed elastic blade are moderate, and of order e (0 10< s< 0.20). Orders of 
magnitude are then assigned to various non-dimensional physical parameters governing 
the rotor blade dynamic problem, in terms of s In the derivation of the governing 
equations, higher order terms (terms of order greater than s) are neglected with respect to 
terms of order 1 term, i.e.. 


0 ( 1 ) + 0 ( 8 ^) - 0 ( 1 ) 

The order of magnitude of various non-dimensional parameters governing this 
problem are given below 

Order 1; 


cos^^,sin<zi^,A^,A, =0(1) 




‘ ®()=^()=0(1) 




d<l> 


‘4ri)=47() = om 


dx^ 


cbcl 


Order 

=0(6*“) 

Order 8 • 
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a rj 


0(s ) 


r r r g' r i 


Order : 


Im„,Im„=0(8“) 


7777’ ^5- 




Order : 


It is important to note that ordering schemes are based on physical understanding 
of the behaviour of actual blade configurations. Hence, care must be exercised in deleting 
higher order terms, based on this ordering scheme. 
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INBOARD SEGNfENT 


a) Schematic model 
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Chapter 3 


Coordinate systems 


The descriptioii of the complex deformation of a rotor requires several coordinate 
systems. The transformation relation between quantities referred in various inertial, non- 
inertial coordinate systems is to be established before deriving the equations of motion of 
the rotor blade. The relation between two orthogonal systems X,, Y,, Zi and Xj,Yj, Zj with 

vectors along the respective axes can be written as 




r - "1 

/V 

•=kl' 




- 


(3.1) 


Where the transformation matrix [T,j] can be obtained using the Euler angles 
required to rotate the j-system so as to make it parallel to i-system. The coordinate 
systems used in deriving the equation of motion for the rotor model are described below: 


3.1 Hub fixed inertial system - R 

The coordinate system ‘R’, shown in Fig. 3.1, has its origin at the center Oh of the 
unperturbed hub The xr axis is pointing towards the helicopter tail and zr is pointing 
upwards. The unit vectors along the three axes are and . 


3.2 Hub fixed moving system - H 

The coordinate system H, shown in Fig. 3.2, is a body fixed system with its origin 
fixed at the center of rotor hub Oh. Prior to perturbational motion of the hub, H-system 
coincides with R-system. If 0x, 0y, and 0z represents the sequential yaw-pitch-roll 
rotations, then the transformation matrix [Thr] can be written as 
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(3.2) 


'l 

0 

0 

COS0^ 

0 

-sin0^ 

COS0^ 

sin0. 

0" 

0 

COS0^ 

sin0. 

0 

1 

0 

-sin0^ 

COS0, 

0 

_0 

-sin0^ 

COS0^ 

_sin0^ 

0 

COS0^, 

0 

0 

1_ 


3 

Since 0^,0^, and O^are assumed to be of order e^, sine and cosine functions can be 
approximated as: 


sin0 «0 ; COS0 » 1 

Substituting the approximation, [2]^] can be simplified as: 



1 

0,0. -0, 
0 , 0 , + 0 . 


0, 

1 

0.0, -0, 



(3.3) 


3.3 Hub fixed rotating system - 1 

The coordinate system 1, shown in Fig. 3.3, rotates about axis with variable 
speed Qq of the rotor. Its origin is fixed at hub center . This system can be obtained 
by rotating H system by an azimuthal angle of the k'’' blade about Zj^ axis. The 
transformation matrix is given as: 



cosi/A* siny/^ 0 
-siny/^ cosy/^ 0 


0 0 1 


(3.4) 


Where, y/^ is azimuthal angle of k‘^ blade 


271 

Wk =V'o+(^-l)Tr 





3.4 Rotating system - 2K 

The coordinate system 2K, shown in Fig. 3.4, is a blade fixed coordinate system, 

which rotates with ^ blade. The origin of the 2K system is at the location of the blade 
root A (Fig. 3.4), which is at a distance +eje^jfrom the hub center. 1 -system and 2k- 

system are parallel. 

3.5 Preconed, rotating system - 3K 

The system 3K, shown in Figs. 3.5. a and 3.5.b, also rotates with blade. This system is 
obtained by rotating 2K-system by an angle (precone angle) about axis. The 

transformation matrix between 2K and 3K systems is given as: 



0 K 

1 0 
0 1 


(3.5) 


3.6 Predrooped, preswept, pitched, blade-fixed rotating system - 4K 


The 4K system, shown in Figs. 3.5.a and 3.5.b, is blade fixed system with its origin 
at pitch beating B (or flex beam-blade-torque tube junction) of the blade. It may be noted 
that the pitch axis of the blade is along direction and the blade reference elastic axis 

is along the direction. While changing the control pitch input of the blade, the elastic 
axis will move on the surface of a cone whose vertex angle is described by the angles 
and as shown in Fig. 3.5.c. 

The 4K system is obtained by the following steps: 

Stepl: Translating the origin of 3K system by a distance ‘e 2 ’ along ^ 3 ^. 

Step 2 : Then rotating the system by an angle (pre-sweep angle) about axis. 

Step3: Then rotating the system by an angle /3^ (pre-droop angle) about y-^^. axis. 

Step4: Then rotating the system by an angle 9 j (pitch input) about . 



The transformation matrix is given as: 


1 



{Ps cosOj+Pj sm0;) 
(-/?^ cose^ + P^ sm0,) 


-{P^ cosOj + Pj sm0/) 
COS0^ 


{Pa cosOj-p^ sm0,) 
sinS^ 
cos 


(3.6) 


3.7 Undeformed element coordinate system - e 

The e-system, shown in Fig. 3.6, has its origin at the inboard node of the finite 
element. The vector , is aligned with the beam element elastic axis; while the vectors 

and are cross sectional coordinate axes. For the straight portion of the blade, the 
i^xe’^ye’^ze) system has the same orientation as system. For the swept-tip 

element, the e-system is oriented by rotating the 4K system about by anhedral angle 

, and then about by the sweep angle . 

The transformation matrix between 4k and e-system is given as: 




^x4k 

^ e 

ye 

II 

A 

^y4k 

Ae. 




(3.7) 


For the element in the straight portion of the blade 



1 0 0 
0 1 0 
0 0 1 


(3.8) 



For the swept-tip element 



cosA^ cosA^ 
sinAj cosA„ 
-sinA, 


-sin A, 
cosAj 
0 


cosAj sinA^ 
sinAj sinA„ 
cosA^ 


(3.9) 


Where, A^ is the blade tip sweep angle, positive for backward sweep, and A^ is 
the blade tip anhedral angle, positive upward. 

3.8 Rotating, blade-fixed system - 5K 

The 5k system, shown in Fig.3.7, Is cross-sectional coordinate system of the 

blade. In the undeformed state of the blade, both e and 5k systems are parallel. But, 
the origin of the 5k system is at a distance from the origin of the e - system. During 

elastic deformation of the blade, the 5k system translates and rotates with the cross- 
section. After deformation, the origin of the 5k system, from the origin of 4k system, is at 
the location given by 


n-1 A 

,<=1 J 




(3.10) 


The transformation matrix between e and 5k system is obtained following a flap - 
lag sequence of rotation. The Euler angles are respectively -/3^. and corresponding to 
the local slope of the deformed blade in flap and lag directions. The transformation 
matrix is given by: 



(3.11) 


■ cos^^ 

sm^, 

o' 

■ cos^. 

0 

sinp,' 

-sin^. 

cos^^ 

0 

0 

1 

0 

0 

0 

1_ 

-sinpk 

0 

cosp,_ 


Since the angle /3i and are of order 0(^s), the transformation matrix can be 
simplified by noting that 

sm0 *0 ; COS0 « 1 

The Euler angles can be expressed in the terms of the local slope of elastic 
deformation of the blade as 




dx 

dx 


n-\ 

Note: x = Y,{le)i+x^ 

Substituting the above relations in the transformation matrix [7^^] yields 



1 

0 




■^k^k 


(3.12) 


3.9 Coordinate system-6K 

The 6K system, shown m Fig. 3.8, represents the cross-sectional coordinate 
system in the deformed configuration of the blade. The term represents the 

directions of the cross-sectional principal axes. 6K system is obtained by rotating 5K 
system about through the angle where ©^represents the geometric twist 

angle of the cross-section and represents the elastic twist. The transformation relation 
is given as 



( 3 . 13 ) 



0 

cos{(f), +0^) 

sin((/>i+0c) 
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Figure 3.4 Rotating hub system - 2 
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b) Top view 
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Figure 3.7 Rotating blade fixed system - 5k 
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Chapter 4 
Kinematics 

During operation, the rotor blade undergoes in in-plane bending (lag), out-of- 
plane bending (flap), torsion and axial modes In addition, the hub center has both 
translational (Rx, Ry, Rz) and rotational (6x, 6y, 0z) motion The formulation of inertia 
operator and aerodynamic operator requires a proper description of kinematics of the 
blade motion In this section, an expression for the absolute velocity at any arbitrary point 
‘p’ on the blade is derived. 

4.1 Position vector of a point 

The position vector of any arbitrary point ‘p’ in the nth finite element of the 
deformed blade with respect to the hub center <9^ , is given by 

n-\ 

+^\^x2k 

1=1 

(4.1) 

Transforming all the unit vectors to the 4K-system and neglecting the higher order terms, 
the position vector can be written as; 

/l-l 

^ = ^x^A-aiJS, cosOj +0^ sin^,) + e, +^2 +^{h), +cosA, cosA^(x^ +u^) 

z=l 

-t-v^ sin A^ cosA^-w^, sin A^ cosA^ cosA^ cos((9q +(f>p)-M>' cosA^ cosA^ 

sin(^G +^fc) +sin A^ cosA^ cos(^o -f-^^)-sin A^sin(^o +^it)] +C[^k cos A^ cos A^ 

sin(^o +^t)~'^k cos Aj cos A^ cos(^q +A)~ sin A^ cos^ sin(^Q + ^j.) - sin A^ 
sin(^G }+ey 4 jiacos&^ +ei(>9,cos^, sin + cos 6>j -i- 

sin ${)- sin A^ {x^ +«*:) + cos A^ + ? 7 [v' sin A^ cos{0q + <2^* ) + si’^ 
sin(^o + + cos A, cos(^g + + Ch^k sin A^ sin(^o +A) + ^k sin A^ cos(^g + A ) 

- cos A^ sin(<9Q + )] } + ~ ^ sin {-P^ cos 6j + sin Oj - p^, cos 9j ) 

+ e^{-Pd cosOj + P^sin$j) + cosK^s\nK^ix^ +“ifc) + Vi sin A^ sin A^ + w^ cosA^ +t] 


) + (?^k + Uk)^xe + ^k^ye + ^k^.e + , 
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[-v^ cosA^ sinA^ cos(0q, +4>k)-'Wk cosA^ sinA^ sin(6(; +^^.) + sin A^ sinA^ 
cos(0c +<l'i)] + ?K cosAj sinA„ sin(0jj +<f>t)-w[ cosA^ sinA„ cos(0c +04)-sinA^ 
smA^sin( 0 g + cosA^ cos( 0 g +<^> 4 )] } 

(4.2) 

In equation (4.2) all the length quantities are non-dimensionalized with respect to 
the length of the blade ‘1’. For the sake of convenience, the non-dimensional length 
quantities are referred without an over bar (~). The expression given in Eq. 4.2 is general, 
in the sense that for position vector of any point in the straight portion of the blade As and 
Aa are set equal to zero. Eq. 4.2 can also be written in symbolic form as 

] (4.3) 

4.2 Angular velocity vector 

The angular velocity vector c5 of k* blade consists of three components. They are; 

• The variable speed (Q) of the rotor. 

• The rigid body angular velocity of the hub due to perturbational 

rotation in roll-pitch-yaw ( 0 x , 0 y , 0 z ). 

• The angular velocity contribution due to the rate of change of control 
pitch input 0.^0; to the blade. 

(£> — + ^0^f^x3k (4-4) 


Where, 


^ rigid ~^0^x^H 


(4.5) 


Where, ( ) indicates differentiation with respect to the non-dimensional time xj /^ , 



Vo=^ot 


Note that 




Transforming the all the unit vectors of Eq. 4.4 to the 4K-system and neglecting the 
higher order terms, the angular velocity of blade can be written in symbolic form as: 


= Qo[a 


+®2^z4iJ 


Where, 



cosBj +6^ sinO; siny/^ -6^ cosyr^ sinO^ +6^ cosdj +0^(-/3^ cos0^ + j8^ sin0^) 


(4.9) 



4.3 Velocity at a point ‘p’ 

The absolute velocity vector V , at point ‘p’ on the deformed beam can be written 
as: 

^ = +rp+i^xr^) (4.10) 

Where, is the rigid body perturbational translation of the hub center , which is 
given as: 


V„ - ] (4. 1 1) 

Transforming all the unit vectors to 4K-system 

Vh = ). \ + (f'« X e,.. ] (4.12) 


Where, 

= C2o^R^(cosif/,^ + cos9j + (3^ sin9^)}+QJRy{siny/^ -cosy^(/3^ cos9j + sin0^)} 

+ QJrXI3p+ Pj cos9j -p^sm9j) 

(4.13) 

i^H)y =^^0^^Jcosv/i(i3, COS0, +i3^ smOp-Py sin0,)-smv/i cos0/} + 
f^(,/^j,{sinv/4(/3, cos0^ + jS^ sinO, - sin0^)+cos0^ cosy/* }+ Qo^^i(siu0,) 

(4.14) 

(F^). = Qo/^^{cosi/r^(- Pj cos9j + P^ sin0^ - p^ cos0^)+sinv^4 sin 6; }+ 

Qo/.R_j,{-cosv/^ sinO^ + sinxf/ p ^ cos9j + P^ sinO^ - cos0^)}+Qo/.^,{cos0^} 

(4.15) 


and: 

rp =e^^JQo[-ae,{-P^ sinO, +p^ cos0^)+cosA, cosA^ +sinA, cosA^v^ 
-cti^ sinA^ +n\-v[ cosA^ cosA^ cos(0^ +Oi)+Vj0^ cosA^ cosA„ sin(0(; +<j>k) 

cosAj cosA„ sin(0(j +<j)i^)~a}[ cosA^ cosA^i^^ cos(0c sinA^ cosA„ 



sin (0c +^i )-4 sin cos(0c }+C{vt cos A, cosA^ sin(0c 

+ cosAj cosA^ cos(0c 4-^^)-©^ cos A^ cosA^ cos { 6 q + 

cos A, cosA„ sin(0c sin A, cosA„ cos(0g^(/»J+ 4 sin A„ sin(0g +<^J } ] 

+ e_^44/Qo[-a0^ sin0; +e^6j{- j3j sin0^ + j8^ cos0^ - cos0^)+e20;(- j3, sin0^ + /3^ cos0^ ) 

-Wa sin a, +Vi cos A, +n{v; sin A, cos(0c +(l>i)-v'J^ sin A, sin(0c +<^>i)+(o; sin A, sin(0c +((>^) 

+ sm A, cos( 0 e + ) - 4 cos A^ sm( 0 c + 0 * ) } + <1 {- v* sin A^ sm (05 + (^ J - v ;(^4 sin A^ cos( 0 c + ) 

+©; sinA, cos(0c +(^i) -©^4 sinA, sin(0c +<^i) -4 cosA, cos(0c +(^ 4 ) } ]+ 444 /J^o [-«4 cos0^ + 
e^OjdSj sin0; + J3j cos0; + J3^ sm0;) + e24(i^rf sinO^ + /3j cos0;) + M 4 cosA^ sinAa+v^. smA^ sinA^ 

+ ©4 cosA^ +«{-V4 cosAj sinA^ cos(0c +(f>t) + v[^^ cosA^ sinA^ sin(0c H-^^)-©^ cosA^ sinA^ 
sin( 0 c +(/»4) -©;4 cosA, sinA^ cos( 0 c +04 ) -4 sinA, sinA„ sin( 0 c +^4) +4 cosA„ cos( 0 c +^4) }+ 
^{4 cosAj sinA^ sin( 0 (j +^4) +V44 cosA^ sinA„ cos( 0 c +^4) -©4 cosA^ sinA^ cos( 0 c +^4) + ©^4 
cosA^ sin A„ sin(0c +04)-4 sinAj sin A^ cos(0c +04) “4 cosA„ sin(0(; +04) } ] 

(4.16) 


and 


coxr = 


'x4k 


CO, 


CO^ 

'■y 


’^z4k 

CO, 


(4.17) 


(Sxr) = e^,, (co^r^-r^(0^)-e^^,(a>^r^-r^coj + i^4,(co^r^-r^o)^) 

(4.18) 


Where, Tx, ry, and ©x, ©y, ©z are the x, y, and z components of 4 and © respectively. 

Substituting various quantities in Eq. 4.10 from Eqs. 4.2, 4.7-4.9 and 4.13-4.18, the 
velocity at point ‘P’ is obtained. This can be written in symbolic form as: 




(4.19) 



The components of the velocity vector are given below 


K = (Fh )x+i>'p)x+ ) 

K = )y + (fp ),, - {C0^i\ - ) 

K- = + (^p ).. + - 0)y>^:c ) 


(4 20) 


(K ) = ^o 4 ^r (cos 5 ^/-^ + sin y/^ {J5, cos dj sin dj )) + (- cos y/^ cos Oj + p^ sin 9j ) + sin 5 ^^ ) 

+ RPPp + P^ COS0J - P^ sindj)- adj{p^ cosOj - P^ sin(9p) + cosA^ cosA^ilj. +sin A^ cosA^v^ 
sm A^ - \p cos{6q + ) - ^ sin(^Q + (j )^ )] cos A^ cos A^ + cos A^ cos A^ + (p^ sin A^ ) 

+ [;; sin(^G + <4 ) + ^ cos(^g + <4 )] ( v;4 cos A, cos A^ - d)[ cos A^ cos A„ - 4 sin A, cos A^ ) - n J 
~ (^) + A cos ] (e, + J (xjt + ) [cos sin A, + sin Bj cos A, sin A J + J 

v^. [sin^p sin A^ sin A^ -cos^p cos A^ j+I^j®;. sin^p cos A„ +[ 7 cos(^g +^p)-^sin(^G + ^j..)](-Vp. 

- j [sin $j cos A^ sin A^ + sin A^ cos 0p ] + j [sin Bj sin A^ sin A^ - cos 6>p cos A^ ]) +[ 7 sin(^Q + <4- ) 


+i^cos(^Q ^^j[sin^P cosA^ sin A^ +cos^p sinA^]+|^jsin6>p cosA^ sin^,. 

\p^ sm 2Bj - pj cos 2(9p -pp\-B^e^ sin y/^ \p^ sin 26 j - p^ cos 2(9p ] - (9^ sin (x^ + )[cos Bj cos A^ 
sinA^ -sin^P sinAj ]-^^Vj. sin^4/'j.[cos^P sinA^ sinA^^ +sin^P cosA^J-^^cW;;. sin^^^P cos^p cosA^ +B^ 
[ 7 COs(^g +^jt)“4^sin(^o +^)](Vfc sin^^(cos^p cosA^ sin A^ -sin^p sin A^)- sin [cos^p sin A^ sin A^ 

+ sin Bp cos A^ ]) + ^^ [7 sm(^(j + 4 ) + ^ cosC^o + A )] (®jfc sin Wk [cos B^ cos A^ sin A„ - sin ^p sin A^ ] - 
sin y/p. cos B^ cos A^ ) + Bye^ cos y/^^ \p^ sin 29 j - p^ cos 2B, - p^ ]+ ByC^ sin y/y^ \p^ sin 2Bi - P^ cos 2J9p ] 
+ By cosy/i^(xi^ +Mj.)[cos^p cosAj sinA„ -sin^p sinA^]+^j, cos^i/'t.Vp.[cos^p sinA^ sinA^ +sin^p cosA^] 
+^y cos^^p.<yj.(cos^p cosA„) + ^y[ 7 cos(^o +4)-#sin(^o +4)](“Cos^^^X(cos^p cosA^sinA^ -sin^p 
sin A J + cos yr^. (cos Bj sin A^ sin A„ + sin Bj cos A ^ ) + [7 sin(0o + 4 ) + ^ cos(^o + 4 )] (~<®fc cosyr^ 
(cosBj cosA^sinA^ -sin^pSinA^) + cos 5 ^/'j.cos^p cosA^ -(^1 +^ 2 )^z(A sin2^P + p^cos2Bj) 

+ B^(Xk +i/^.)[sm^P cosA^ sinA^ +cos^p sinAj]+^ 2 Vj.(sin^P sinA^ sinA^ -cos^p cosA^)+^^ 0 ;i. cosA^ 

sin ^P + [7 cos(^G + 4 ) “ ^ sin(^G + 4 )] (“^ (sin Bj cos \ sin A^ + sin A^ cos ^p ) + sin Bj sin A^ sin A^ 
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- cos dj cos [n sin(^o + ) + ^ cos(^q + (p^ )\{-o>[ (sin 9^ cos sin A^ + sin A^ cos 9,) + sin 9j 

cos A^ 2^/ “ Pd cos 29 j ] + 9^ (x^ + )\p^ (cos 9j cos A^ sin A^ + sin sm A^ ) + (sin 9j 

cosA^ sinA^ -cos^^ sinA^) ] + ^^v^[/?^(cos^^ sinA^ sinA^ -sin^^ cosA^) + PPsm9[ sinA^ sinA^ + 
cx3s9j cosAj]+)9;iyjt cosA^ [^^cos(9; +y5jSin(9;] + ^; [;7 cos((9g +^4)-^sin(6(j +^zij.)](-v; [ppcos9j 
cos Aj sin A^ + sin 9j sin (sin 9j cos A^ sin A^ - cos^^ sin A^ ) \+[p, (cos9j sin A^ sin A^ - sin 9j 

cos A^ + p^ (sin 9j sin A^ sin A^ + cos 9j cos A^ )] j + 9j [;; sin((9Q +(p,^) + ^ cos{9q + <4 )] (-o^ [ P^ (cos 9j 

cos A^ sin A^ + sin 9^ sin A J (sin 9j cos A^ sin A^ - cos 9j sin A^ )] + cos A^ (y?^ cos 9j+ P^%m.9j)'j^ 

(4.21) 
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i^y ) “ ^o4A ; ifis (^0^9 j + sin^i - y?p sin^; ) - siny/^ cos9j )+Ay (cos^/^. cos^^ + siny^j^ {/3^ cosOj + 
Pd sin^j - pp sin^j) ] + i?^sin^^ -aOj sin^j -^^x^p-Ps sin^; + P^ oo%9j - P^ cos<9j) + e 2 ^/(“;^jSin^; 
+ Pd cos 6,)- Uf. sm A^ + Vj. cos + [ 7 cos(9q + (f>P)-^ sin(^Q + sin A^ + sin A^^^^ ) + 

+(!>,)+ ^ cos(^G + 4 )] sin A^ + sin A, - 4 cos A^ ) + i 7 sin 6 »; - + (^) (e, + ^3 + 

S"' )cos (93 -(^)(:>Ci: +Mi)cosAjsinA^(/?^ +y 5 ^cos (93 sin 6 ' 3 )-cos (93 cosA^ ]-(-^)w^ sin A^ 

2=1 

[sinA^(y?p + p^ cos9j -P^sm9j)-cos9j cosA^J-(-^)tyj,[cosA„(y9p + p^ cosBj - P^sin9j) + s,inA^ 
cos 9, ] - [7 cos((9o + 4) - ^ sin(6>o + ^t)]{-v' cos A^ [sin A^ (P^, + p^ cos 9- - p^ sin 9j)- cos 9j 
cos A J + sin A^ [sin A^ (P^ + cos sin 6*; ) - cos cos A J} - [7 sin(^Q + # cos(^g + 

<^ir)]{-®itCOsA^[sinA^(^^+^^cos6'3 -y5^sin(93)-cos6>3COsAj+[cosA^(>5^ +y9^cos6>3 -y^^sin^^^) 
+ cos^3 sinA^ ] ]-9/isir\.9j{—cosy/^ + p^ cos9j sini/Tj^ + p^^ sin^^ sin^2/’^)-^^ej[cos{2/'j.(-/7^ cos^^ + 


p^ sin 9j - Pp cos9j ) - sin yr^ sin 9j ] - e^9^ [cos^^^y. {-P^ cos9j + P^sm9j)- sin yr,^ sin 9j]+9^'^u^ sin y/’;^ 

2=1 

sin ^3 -9pxy. +Mjj.)[cosA^ cosyr^ sinA^ -sin^j sin^^^^t cosA^ cos [sin A^ cosy/^. sinA^ -sinA^ 
sin (9; sin cos A^ ] - 9^(0,^ [cos y/^ cos A^ + sin 9j sin y/^ sin A^ ] - 4 [27 cos(^q + ) - - ^ sin(^c +<t>k)\ 

cos Aj (cos y/^ sin A^ - sin 9^ sin y/^ cos A^ ) + sin A^ (cos y/^. sin A^ - sin 9^ sin y/^ cos A^ )} - 9^ [7 
sin (^(3 + 4 ) + cos( 4 ; + 4 ) ]{“ cos A^ [cos^ 2 /^jt ~ siii ¥k cos A^ ] + {cos^k cos A^ - sin 6^ 
smy/^ sin A^ ) }+ 9yasvn.9j {siny/^ + P^ cos9j cosyr^ + p^ cosy/^ sin^^ ) + 4^1 {Pd cos^j sin^/jj, - p^ sin^j 
sin + Pp cos 9j sin - cos y/,^s\r\9j)->p 9 y {p^ cos 9j siny/^- P^ sin 9j sin y/^ - cos y/^s\n9j)- 9 y 


«-l 


cos^i^i sin4 sin A^ +sin^3 cxis^k cosA^)-4vy,(sin A^Csm^i/y^ sin A^ + 

2=1 

sin 9j cos y/^ cos Ap)-9yCO^. (sin y/^ cos A^ - sin 9j cos y/,^ sin A^ ) - 9y [7 cos(^q + 4 ) “ ^ sin(^Q +A)] 
(-v^ cos A^ (sin y/^ sin A^ + sin 9 j cos y/^. cos A^ ) + sin A^ (sin yf^ sin A^ + sin 9 j cos 52^^. cos A„ )] - 4 [27 
sin(^G + 4 ) + 4 " cos(6>g + 4 ) K- m'y. cosA^(sin^2^jj, sinA^ +sin^3 cos^2^^ cosA^) + (sin^2/^j. cos A^ - sin 4 

cos^i^^sinAJ }- 9 ^a{p^-ppSmei) + d^ e^+e^+Yp*, cosGj -epX;^+uP)cosAj^mKppp+ P^ 

\ /=i y 

cqs 9 j - y^^sin^;)-cos0j cosA^) ]-9^v^^ sinA^[sin A^(y9p + P^ cos9j - P^sm9j)-cos9j cosA^]- 
9^0),^ [cos A^ cos sin ) + sin A^ cos 29^ )] - 4 [7 cos((9o + 4 ) “ ^ sin(^o + 4 )] {-v^ 

cosA^(sin A^(;ffp + p^ cos9j -y9^sin^^)-cos^j cosA^)+sin A^(sin A^(;ff^ + P^ cos9j - P^ sin 9j) 
- cos 9j cos A„ )} - 4 [27 sin(43 + 4) + 4^ cos(4j + 4)] {-<4 cos A^ [sin A^ (pp + p^ cos 4 - p^ sin 4 ) - 
cos4 cos A^ ]+ (cos4 sin A^ + cos A ^ ifip + p^ cos4 “ A sin 4 )) }+ 4nfsin 4 + ^i^xPp cos4 + 


Y^‘Pii-Pd cos 4 +p, sin 9j)-Gj (xj, + ) cos A, [cos A^ {P^ cos4 - >9^ sin 4 ) + sin A^ ] - v,.4 

1=1 
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sm [sin + cos A^ 03, cos 0^ - j3^ sin 0j )]-Q),0j [cos A^ + sin A^ (J3^ sin 0j - p, cos0j )]- 0j 
[^ cos(( 9(3 + <4 ) ~ sin(0Q + (j)!^ )] [-v^ cos A^ [sin A^ + cos A^ {p^ cos 6^ - P^ sin 0j )] + sin A^ [sin A^ + 
cosKpp^ COS0; - P^sin0j)]^-0i [7sin((9Q + 4) + ^cos((9q cos A^ (sin A^ - cos A^ 

(P, sm 0 j- Pj cos 6 j )) + (cos A^ + sin A^ {P^ sin 0 j - p^ cos 


(4.22) 

^z)=^\K^^Wk{~fid + Ps sin^/ -Pp cos6j)+smy/^ sin^J+.^[-cos{i/-^, sin^; +smYk{-pa cos$j 

+ P^ sin 9j - p^ COS0J ) ]+ ^. COS0J - a0; cos0j + e^Oj (P^ sin 0j + p^ cosOj + p^ sin 0^) + e^0i {p^ 
sin(9j -\-p^ cos(9;)+/4 cosA^ sin A^ +Vj. sin A^ sin A^ + aij. cosA^ +[;7cos((9o +(j>P)-^%m(0Q +^P)\ 

{-v[ cos A^ sin A^ - cos \ sin A^ + ^ cos A^ ) + [7 sin((9G +<j>P} + ^ cos(^o + j^)] (v ' 4 cos \ sm \ 

- cosA, sin A^ - 4 sin A, sin A, ) + [^{p, + P^ cos^, )- +e^+ ^u, 

[sm A j (Pp + Pj cos 0j - P^ sm ) + cos A ^ cos A„ sin ^ J + [cos A ^ (p^ + p^ cos 9j - P^ sin 0j)- 

sinA, cosA, sin 4 ] + (^)tytsin A,sm^, +(^)[ 7 COs(^o + 4 )- 4 ^sin(^o + 4 )]{v; [sinA^^p + cos0j 

- p^ sm 0j) + sin 0j cos A ^ cos A^ ] + cos A^ (/?p + p^ cos 0p - p^ sin 0j)- sin A^ cos A^ sin 0j } + 

[V sm(0o + 4 ) + ^ cos(4j + 4 )] [sm A^ (>9p + p^ cos 0j - p, sin 4 ) + cos A, cos A^ sin 0j'] + sin 0, 

sm A^ }+4 :^cos4 ('^Wk ~ Ps ^^^¥k cos4 -Pd ^^^¥k sin 4 ) + 4 :^i(sm?i/'j. cos0i + P^ cosy/,^ cos0j + 


sin^/ -(§K +^k) 


n~l 

P^ cosy/,^ sin 0j - p^ cosy/^. sin 4 ) + 4 ^ 2 (A cos^i^^^j, cos^ + Pd c.osy/,^ sin 0^ + cos0j sin y/^) + 0P^'u, 

i=l 

sm^f^^.cos4 +dpx^. +M^.)[sin^i/'^. cos0j cosA^ cosA^ -cos^^/’j.smAj]+4vi[cos^i/'4. cosA^ +sin^ji. cos^ 
sm A^ cos A^ ] - 4<y^ sin y/^ cos 0j sin A, + 4 [?7 cos(^g + 4 ) “ ^ sin(6(j + 4 )] {v;(cos y/^ sin A^ - sin y/^ 
sm y/^ cos A^ cos A^ ) + (cos y/^ cos A^ + sin y/^ cos 4 sin A^ cos A^ )} + 4 [7 sin(4 + 4 ) + “s cos(4 + 4)] 
{a)[(cosy/^ sinA^ -sinyr^ cos0j cosA^ cosA^) -sin^/j, cos^; sinA„}+4«cos^;(sin54/'j. + p^ cosy/,^ cos4 

Pj sin 4 cos Yk) + ^Ps ¥k cos 4 + Pd sin ¥k sin 4 ~ Pp sin Yk sin 4 - cos cos 4 ) + 4^2 


n-l 

\P, sin 5 /^ cos4 + Pd sin^^*: sin4 -cos^/^, cos4]-4v2"' ^^Wk cos4 -Gy(x^ +M*)[sin 5 i^i, sin A^ + 

1=1 

cosYk cos4 cosA^ cosA^ ]+4vfc[sinv^^ cosA ^ -cosYk cos4 sinA^ cosA„]+4 ^^Wk^k cos4 sinA„ 

+4 [0 cos(4; + 4 ) “ # sin(4 + 4 )] {^k (sin ¥k sin K + cos Yk cos 4 cos A^ cos A^ ) + (cos A^ sin Yk ” 
sin A^ cos A^ )} + 4 [?7 sin(4; + 4) + ^ cos(4 + 4 )] (sin Yk sin \ + cos Yk cos 4 cos A^ cos A^ ) + 

(cos?^/'^ cos4 sin AJ ]+0^a{p^ +4 cos4)-4(e, +^2 -^si^k +«;t)[sin A^(4 +p^ 

1=1 

cos4 -4^sin4) + cosA^cosA^sin4 ]+4vi[cosA^(yffp +y^^cos4 -j4sin4)-sinA^cosA„sin4] 
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+0,6),^ Sin sin 0j + 6^ [;; cos(6>g + ) - ^ sm((9o + ^^)\ {v' [sin cos 6j - p^ sm Oj ) + sin 9j 

cos A^ cos A^ ] + cos A^ {p^ + p^ cos 0, - p^ sin 0j )- sin A^ cos A^ sin 0j] + 0. [;/ sin(( 9(3 + ^ ) + ^ 
cos(^Q + )] [sin A ^ (Pp + Pj cos 0j - p^ sin 0j) + cos A^ cos A^ sin 0/^ + sm sin A^ | + 4^ cos 0; - 

n-\ 

4 iPp sin 0j)-0j iPs COS0J + p^ sin 4 ) “ 4 + ‘h )[sin + cosA ^ cos A„ (P^ cos0j + P^ sin 4 )] 

z=l 

+4v;i.[^sAj -sinA^cosA^(y4cos4 +y4sin4)] + 4‘^fc(si*^A^(j4cos4 + j4sin4))‘''4 [7*^(4? + 4) 

sin(45 + 4 )] [ sin A^ + (4 cos 4 + >4 sin 4 ) cos A^ cos A^ ] + [ cos A^ - (4 cos 4 + y4 sin 4 ) sin A^ 

cos A^ ]^ + 4 [7 sin(4j + 4 ) + cos(4; + 4 )] [sin A^ + cos \ cos A^ (4 cos 4 + >4 sin 4 )] + [sin A^ 

(4 cos 4 + 4 sin 4 )]}] 
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Chapter 5 

Equation of Motion for Rotor Blade 

The coupled equations of motion of the rotor blade can be derived using 
Hamilton’s principle The mathematical form of Hamilton’s principle is stated as follows; 


\(5U-6T-5W^Yt = Q (5.1) 

Where U is the strain energy; T is kinetic energy; is the work done by the non- 
conservative external loads 

In this chapter, the expressions for the variation of kinetic energy and strain energy of the 
rotor blade are derived 

5.1 Kinetic Energy of Blade 

The total energy of the beam, T, is defined as; 

^ * 0 

Where V is the velocity of an arbitrary point ‘p’ on the blade cross section with respect 
to the inertial reference frame. The variation of kinetic energy, ST can be written as; 


ST = £ pV.SVdr]dCdx (5.3) 


Substituting for the velocity V from Eq. 4.17 and integrating ST by parts with respect to 
time, yields 
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(5.4) 


i'A 

^T = Yu I \\,p[^uSih + ^M+2.5w,+Zl5V,+Z'Jw[+Z^6(l>,']iT]d<;dx 

^ 0 

Where the terms Z„,Z,,Z„,Z',Z; and Z^are the coefficients of 5ii„6v„dw^,dv[,5w[ 
and in the variation of kinetic energy expression. 

Integration of the expression over the cross-section yields: 


(O, 

6T = J [Z.*, +Z>. +ZJw, +Z’M+z:Sw', +2,5^,}* (5.5) 

i 0 

After eliminating the higher order terms, the expressions for Z„,Z„,Z„„Z',Z^ and Z^are 
obtained. 

The cross-section integrals are defined as: 

m = \l^pdrid^ 

m„ = \\^pndndi; 

Where m is mass per unit length of the blade; mr]^ and are the first moments 
of cross-sectional mass per unit length; are the cross-sectional mass 

moments of inertia per unit length of the beam. 

5.2 Strain Energy of Blade 

The formulation of strain energy outlined in this section essentially follows the 
procedure given in the Refs. [6], and [5] 
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5.2.1 Strain Energy of the beam Element 


The Strain energy of the beam element is given by 


u 




" •> 

T 




O’ 

XX 


XX 

Yxn 

> < 






dr] di^ dx 


(5.6) 


5.2.2 Explicit Strain-Displacement Relations 

The expression for non zero strain components written in terms of u, v, w and ^ 
are given in Refs. [6], and [5]. They are: 

+«xV 

+ aT„ ) - [tj cos (e^+(l>)-C sin (0^ +(/>)] 

- T] sin (0g + 0 ) + C cos (e^ +(j))]w^+ri (f^ ^ ) + C + T^ofxr , ) 

Yxn =f:cn +CC¥r^-C{(l>x+(l>o) 

Yxc =Yxi; +avc -Tl{(l>x+<f>o) 

Where, 

t={vxx cos sin 0c ) (-v^ sin 0^ + w, cos 0^ ) 


(5.7) 



The underlined term in represents the axial strain at the elastic axis. These 
strain expressions can be simplified using the following assumptions; 

• The transverse shear at the elastic axis is assumed to be zero. 

• The warping amplitude a is assumed to be equal to -(p ^ . 

The simplified strain components can be written as 

+C^)<l>x -WK - [^0 (^V^n 

- [t 7 cos (0g + 0;, ) - C sin (0^; + )] - [n sin (0^ + ) + C cos (0^ + (p^ )] 

r:cn=-{¥rj+^)t-^<l>o 

Yxi =-(v^f -r])(p,+rit 
where 

t = i^xx cos 0G + sin Og ) (-v^ sin 0^ + cos 0^ ) 

(5.8) 


5.2.3 Stress-Strain Relations 


Assuming that the blade is made of isotropic material, the stress-strain relationship is 
given by the following equations: 


r ■> r- -n r ^ 


^ XX 


E 

0 

0 


^ XX 

^ XT} 

> = 

0 

G 

0 

< 

Txn 



0 

0 

G 


Yxi J 


(5.9) 



5.2.4 Strain energy variation 


The variation of strain energy of the i*’’ element is given by 





T 

f N 

'^Yxn 

> < 

XT] 

v ^ J 


^ 

b 


dr] dx 


(5.10) 


The variation of the strain components are given as follows: 

5 + W/ + (n ' + C ' ) ^Mx - -[z,{CWrj-mz)] S4>x 

- [l cos (ec+(t>K)-^ sin (0c + (l>^ )] (5 +^5w^+ wj(l> ) 

- [ri sin (0c + (^^ ) + C cos (6^ + )] ) 

5r^=-(v,+^)s^x-^st 

Yxi =-(v^c -'n)d(l)^+r]5(l)^ 
where 

d(j)^ = (-5 sin 0g+5w^ cos 0 g ) (v„ cos 0c + w^ sin 0c ) 

(-v^ sin 0c + cos 6^ ) (.5 cos 0^ + 5 w^ sin 6^ ) 

(5.11) 

It is assumed that the variations of strain components are of the same order as the 
corresponding strain components. Substituting the above expressions in the strain 
energy variation, we get the components of variation of strain energy as follows: 


«^=l{ 

5m, +m, +0.5v,' +0.5w^ + 0.5(77' +C)<Px -[^ 0 (CWr, -l¥^)](l>x 

-[t] cos (0c + ) - C sin (0c + )] V;„ “ [n sin {9g+<I>k)+^ cos (0c + (p^ )] Wxx ]} 



5v^ {-rjGsinO^ (v^ cos0^ + sine^)^-(j>^ (-tj ) 

+(v^ COS0C + sin0c) (-v^ sin0^ + cos0p)r] ] 

+C G sin (v^ cos 0^ + sin ) [-(/.^ (v^n + C ) 

-(v„ COS0C + sin0g) (-v^ sine^ + cosO^ )C ] 

+ +“x + 0-5v^ +0.5w^ + 0 . 5 ( 7 ]^ + “[^0 (CVn -nv'f )] 0 ;c 

-[t]cos( 0^ +<^>^)-C sin(0c +<j>^)'] -[i7sin(0^ +<I>k) + C cos{e^ +(I>k)]^^]} 

+ 

5v^ {nGcose^ (w^ COS0C -v^ sin0c)[-<^;, (-?? +V/^ ) 

+{w^ cosGc - sin0o ) (v^ cobSq + sinS^; )] 

+CGcos0e (w, COS0C -v^ sine^)[-0^ (y/^ +(^) 

-(w, COS0C -v^ sin0e) (v^ cosO^ + sinO^)] 

+E (-77 cos(0Q +(p)+t^ sin(0^ + ^)[- 0 „V 7 + + 0.5(v^ + ) 

-<^;cTo(V'„C -V('n) + 0.5(t)Xri^ +0 

- (C cos(0e + <^) + r? sin(0c +</»)) 

- v^(77 cos(0c +<f>)-^ sin(0c +<!>))]} 

+ 

5 w, {77 G cos 0G ( cos sin 0^ ) [-(^^ (-T 7 + ) 

+ 77 (w^ COS0C - sin0c) (v„ cos0c + sin0c )] 

-CGcos0c (v^ COS0G + sin0c)[-(/»^ K +C) 

-C (w, cos0^ - sin0c) (v„ cos0g + sin0c)] 


+Ew^ 


-VK +«;. +0-5(vf +w") + 0-5(r7' +C ')</>' -[to(Cv^^ 


- cos(0G +0) +77 sin(0e +(j))) 

- (77 cos(0^ +</►)- C sin(0c + (^)) ]} 



{77Gsin0c {w^ cosO^ -v, sin0^)[-<^^ (-77 +xj /^ ) 

+71 (w, COS0G -V, sine^) (v^ cosO^ + sinO,,)] 

-C G cos 0 ^ ( cos 0 ^ - V, sin 0 ^ ) [-(j)^ + C ) 

(w^ COS0C - sin0^) (v^ cos0e + sin0c )] 

(C cos (0^ + (^ ) + T] sin (0G + 0)) + 0.5 + wj ) 

-4>.^o(¥n^ -Vi;Tl) + 0.5(j)Xri^ +C^) 

- cos( 0 c +(^)+77 sin( 0 ^ + 0 )) 

-v^(T 7 cos( 0 e +(j))-^ sin(0(, + 0 ))]} 

+ 

^<t>{^[vxx (C cos( 0 e +(^)+T]sin( 0 (; +^))-w^ (i?cos( 0 (j +^)-C sin( 0 (- +(^))] 
[-0^1// + + 0.5 +wl)- (l)^x^ - 17/^17 ) + 0.5(^^' (t7 ' + C ' ) 

-w^ (C cos(0o +(l>)+ri sin(0o +<^))-v^ (rj 005 ( 0 ^ +<f>)-^ sin(0^ +'/'))]} 

+ 

^(l>x {G(ri -y/^ )[-(l)^ (-T7 +V'J+T7 COS0C -v, sin0c)(v^ cos0^ + w„sin0(j )] 

+G (-V/, - C ) [-</';c K + C ) - C COS0C - sin0(; ) (v„ cose^ + w^ sin0(;)J 

+E (-To (v/^C ^ ' ))[-^xx¥ +u,+ 0 . 5 (vl+ ) 

-'^o<Px (Vr,^ -Wil) + 0-5<^ ' ( 1 ? ' + C ' ) 

(c cos (0(; + (^ )+ T7 sin (0^ +(/.))- (t 7 cos (0c + (|> ) - c sin (0^ +<!>))]} 

+ 

5(j)^ ^-Ey/ [-(^^ 17 / + + 0.5 (v^" + ) -To(^^ (i/a^C -¥iV) + 0-5</»" (r?' + C ' ) 

(C cos (0c + (^ ) + 77 sin (0c + (^ )) - (t 7 cos (0c + (^> ) - C sin (0c + <^ ))]} 

^dT]di!^dx 


( 5 . 12 ) 


Now we define some cross-sectional integrals associated with the strain energy: 



EA= \\[E)lr]di; 

EAn. = jl[Er,}lrid^ 

EJC, = fj[E(}lndC 
EI,^ = jlEC-}lr^d( 

Ela = JJ[£t)"]rft,dC 

^'^,c = ll[Eni;]i’idC 

EAC, = /][£(>)“ +C’)]rfTjJC 
£^C, = Jj[£r,(,7'+f=)}i,,<iC 
£^C, = |J[£C(n=+C»)]rf,,<iC 

E^^o = jjlEvjdrtdi; 

EADj = jj[Eljij/^r)d^ 

EAD, = fj[EC^p,jdi 
BAD, = jj[E,y^]l,^d( 

EAD, = j^E,l,(,^^+C)}l’ld( 

EAD, = JJ[£v< {CVn-W; )'}lndi 



EAD^ = J|[£'(CV', -VV; 

£iD; = jl[£C(cv',-i?n)]it)df 

EAD[ = j}[£(7,“ +C = )(Cn 

G/„ = jj[G{v‘+i^)yr,d( 

GJ, = jlG(^W,-riVi)]lndi 
GJ, = jj[G{v;+Yi)yvdi 

Certain additional terms are defined as : 


EAt]^ = EAq^ cosO^; -EAi;^ sinQ^. 

EAd;^ = EAt]^ sinOc + EA^^ cosO^ 

EAQ = EAQ cosOjj -EAC^ sinfi^ 

EAC2 = EAQ sinGg + EAC2 cosB^ 

EAD^ = EADy cos0(j -£L4Z), sinO^ 

EAD2 = EAD^ sin 6 (j + EAD2 cosOq 

EAD[ = EAD[ cos0£j -EAD '2 sin^^, 

EAD 2 = EAD[ sin^c + £^Z)' cos0g 

Elr^Z = COS0G 51110^ 

= Elj.^ COS0O sin0G 

EIi-^ = EIi-i- sin0(; cosQ^ 

GJ = GJ^+ 2 GJ,+GJ 2 

After the integration over the area, we get the following line integral after separation of linear and nonlinear 


terms: 



linear 


5u^ {eAu^ -( j)^EADa -(j>,ToEAD'^-w^EAi;^ -Vj^EAq^ } 

+ 

5 m , \q.5EA [v] +wl) + 0.5<l>lEAC, -(j> ( w„;^+ v„£^)J 


nonlinear 


+ 

5vjv„ 

-W, 


EA^u^+ 0.5 (^v^+w^'j^-(l)^EADa+ 0.5(1)^ EAQ -^A^^EAD'^ 


[eai;, +(j)EAn ^ ) - (^EAn^ -<I>eai;^ ) 

+sin9c (v„ cosec + w,, sin 0^ ) [0, (-GJ^ -GJ.)-</»oGJo]} 


nonlinear 


+ 


{KEAD, -u^EAt]^ +(j)^T:,EADl + w,, (£7^^ cosO^ -EI^^ sinO^) 

+^xx{EIi^ COS 0 C -EI^^ 


' linear 


5v„{cos0<;(iv,cosfl5-v,sine<;)[^,(ffi/„+a/|)+^„a/„] 

-0.5(v=+w^)li^;-0.5,J=Sq+^[-^„IiA+(», + 0.5(v,"+«^))SiC 

- w, [(i^ + S^7) sin a„ + (1 ^ - hJ) cos e„ ] 




+ 


{EIi;t; EI^^ j sin 0 £; + [eI^^ + EI^/. j cos 0 g | 


nonlinear 


6w, w. 


EA[u,+0.S(^i‘, + Mf))-,j,^EAD, + Q.5^lEAC„ -i.xfiAD'^ 


[m. +^EAn,)-v^{EAn. -<(.£<.) 

cosOg (v^ COS 0 C +M'^sme(j)[ 0 ^(Gyo +G/i)+^oGyo]j 


nonlinear 


{^xxEAD^ -u^EAI;^ +(t)^T^EAD2 + cose^ +EIi.^ sinBc ) 

(Elr^i; COS 0 O +EI^^ sin 0 e)}^^ 


/ linear 



5w^ {sinOa {w, cos0e -v^ sin0^)[0^ {G^o + GJ,)+(I>oGJ,] 

-0.5 (v^ +w^ )'^^Co - 0.5^1 EACj +<|)^-(j)^^EAD^ +0.5(v^ +w^^^EAri^ 

-(j)^T, EAD; + 0.5(t>^EAQ - + El^'j cosO^ + (-£7^^ + sine^ ] 

-Vx.[(^^cc -(£^ + ^)sin0c]} 

nonlinear 

+ 

S(t> {v^ + («. + 0-5 (v,' + ))£^ -<i>,ZoEAD^_ + 0.5 f;EAC^ 

~^xx l^'^Cn '^^nn ^®®^g J~^xc sin 0 g + El^^ cos 0 g J 

+(f> [-(^«^ + (m, + 0.5 (v,^ + ))^] -<I>,t,E^[ + 0.5 EiQ 

+^^C^)cos0c +(-£/^^ +£J^^)sin0c] 

-V;„[(^ - ^) cose^ -{Ei;^+Ei;^) sin0^] 

+^A-<I>xxEAD^ + (m^ + 0.5(v^ + wl)EAr\^ -(I>^t^EAD[ + 0.5EAQ(j); 

+^x. cose a -EI^^ sin0(,]+v^[£/j^ cos0^ -£^sin0c] 

+0[-(/»^£L4i)2 +(M + 0.5(v + w))£^Ca -<t>TEAD^+ 0.5 f;EAC^ 

-wJi^EI^^ +EI^i.)cosec + [-EI^^ +£/^^)sm0c] 
r -EIr,r,)cosec -( eI^^ +EI^i. )sin0(;l} 

L- ^ ^ ' -sj nonlinear 

+ 

S4,, {-GJ^,^ -V, [-4,^EAD, +EAjy^, -4,,T^A^-w^^,-v^E^l\ 

L- -J/ Imeai 

+ 

{-To [o.5£^£; (v,^ +wl) + 0.5EAD'£ -0 [w^EA5[-vjADi)_ 

+0, [-<I>^EAD, + [u^+0.5[vI + wI))eAC, -< l>^T,EAD:+0.5(j>^EAC, 

-w^ (eAC^ +(f>EACA - (eAQ + <1>EAcM 

' ' ' nonlinear 

+ 

50 „ {^„EAD, - EAD,u^ +(t>.roEAD, + w^EAD, + 



50XV \-^-^EADq (yl + wl 'j-O.SEAD^^^ +^[-w^EAD^ -v^EAD, )} 

^ ^ J nonlinear 

^dx 

( 5 . 13 ) 

The underlined terms are those, which can be Imearized by replacing the extensional strain term 
by the axial inertial force. This will be shown m the later part of the report. 



Chapter 6 


Formulation of element matrices associated 
with kinetic energy variation 

6.1 Finite element discretization 

The variational expressions associated with the kinetic and potential energy of the 
rotor blades is nonhnear. The unknowns are the deformation functions u, v, w and (p. 
These are dependent on both space and time. The spatial dependence is eliminated using 
a Rayleigh Ritz finite element formulation. The blade is divided into sub-regions (finite 
elements), shown in Fig. 6.1, and the total dynamic potential is calculated for each sub- 
region. By applying Hamilton’s principle to each sub-region, a discretized form of the 
equations of motion can be obtained. In this development, each sub-region is modeled by 
a straight beam type finite element. These beam elements are located along the elastic 
axis (line of shear centers) of the blade. 

The discretized form of Hamilton’s principle is written as: 


]j^(SU,-5T-Swr,)di^0 

k '=1 

Where, 

N = Total number of finite elements in the model 
6U^ = Variation of the strain energy in the z'* element 
dT^ = Variation of the kinetic energy in the z* element 
= Virtual work of the external loads in f' element 
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The degrees of freedom for the beam finite element are discretized in space and time in 
the following manner 

-v1 rw" 0^0 0 

w _ ^ ^ 0 {W} 

^00 {<^^f 0 {O} 

L“J [ 0 0 0 

where, {<z), f , f are space dependent interpolation functions; and 

[V } , {W } , {<!>} , {U } are time dependent nodal values of v, w, <j> ,u respectively. 



The nodal coordinates are shown in Fig. 6.2 

The variation of the displacement functions for the beam can be written as; 

~Svl 0 0 0 Hsvy 

_ 0 0 0 [gw] 

S<f> ~ Q 0 {(f>yy 0 {^} 

L‘^"J [ 0 0 0 {^.f 

In this development, Hermite interpolation polynomials are used for the space 
dependent interpolation functions. A cubic Hermite interpolation polynomial, is 

used for the bending deflections (v,w) and a quadratic interpolation functions, is 

used for the torsional rotation {(f) , and the axial deflection (u). The mathematical 
expressions for these polynomials are given as; 
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W = W = {4>c} 


1 - 3 ^^ + 2 ^^ 

‘.{-e+4‘) 


( 64 ) 


W = {A} = {i>,} = 


l-3^ + 2^^ 
4#-4f 
-^ + 2f 


( 6 . 5 ) 


Where: 



= Span-wise (axial) coordinate of the beam element. 
= Length of beam element. 


For bending deformations, the nodal parameters are the displacements and slopes 
at both ends of the beam element. Therefore, the resulting elements will have inter- 
element continuity for both displacements and slopes. In addition, because of the cubic 
Hermite interpolation polynomial, bending strains vary linearly over the element length 
The quadratic interpolation functions are used for torsional rotation (^) and the 
deflection (u). This polynomial has the capability of modeling a linear variation of strains 
along the element length and therefore provides the same level of accuracy as the beam- 
bending element. The nodal parameters for these elements are chosen as the values of the 
displacements fimction at the two end nodes and at the mid-point of the element. 
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The resulting beam element has 14 degree of freedom: 4 in-plane (lag) bending 
degrees of freedom, 4 out-plane (flap) degrees of freedom, and 3 degrees of freedom each 
of torsion (^), and axial deflection (u). The nodal degrees of freedom are shown in Fig. 
6.2 


6.2 Element matrices associated with kinetic energy variation 

The beam element matrices associated with the kinetic energy variation are 
derived by substituting the assumed expressions for the displacements functions (eq. 6.2) 
in the kinetic energy variation 5T (eq. 5.5) and carrying out the integration over the 
length of the beam element. The resulting variation of the kinetic energy can be written in 
the form: 


X' 


X' 

xJ 


K 


+ 




X' 




L Jl4x3 


X. 


e. 




( 6 . 6 ) 


Where represents the vector of unknown nodal degrees of freedom 



(6.7) 


50 


Detailed expressions for the various matrices defined in Eq. 6 6 are given as follows: 

6,2. 1 Mass Matrix \M \ , ^ : 

T 

-l]xa^, sin(^^ +0 Jcos(^^ +<D^ 

h 

\M\ 2 \ ~ f +^ir)cos(^Q +^jc )} + In^,4-[cos^(^Q +^ji:) 

-sm-(^^+(D^)]}{<D;}{cD:f]^ 

rf r1 

[^13] = -/ [m^n, Sin(^G + ) + ^Cm COS(^G + K )} {^c } {<!> J 

^C~ T~ 

[^14 ] = -J {m^ cosiOa +^1:)- <n, sin(^G + )} {<1^: }{<!',} ^ 

n L 

[^21 ] = J [{( ) cos(^G + ^k) + Ini;,, -[ cos' (Oo + Ojc ) - 

0 

sin=(e„+<I.,)]){<I>:}{<I>;}"]<fc = [M„f 

[M22 ] = I l^w {<D, } {O, f -f {im,;, cos' (Oq -t- ^>j5. ) -I- Jm^ sin' ( 0 ^ -f ) + 

n 

2 sin(0o -i- ) cos(^e + )}{^: }{o; f dx 

!<r j.-| 

[M,3] = J {/w7„cos(^G+^ic)-"*Csin(^G+‘I’ir)}{‘I’c}p9} ^ 

{) L 


[M„] = -j {m;7.sm(6i„+4>,)+m<,cos(6lo+®;t)}{<I>;}{4>,} dx 




. 141811 




[W„] = [0] 

[M„] = [W„f 

K.! = [0] 

,j L J 

6.2.2 Damped Mass Matrix[A/‘'] ^ : 

[Wn] = j w [- sin cos cos A^ - cos A^ sin A^ cos A^ + sin A^ sin A^ cos A J (O J {<I>^ 

cos(<9<5 +<t>k)+^m siii(<9o +(j>k)\ {cosA^ cos A^ +sin A^ sin A^ cos A^}{0 jjo/j dx 
[^2 ] = ni cos A, cos A„ j"” dx 

['^•0 = (4. + + A)oos(«<, + 

- cos A., cos A (cos- {Oo + <^, ) “ sin^ (^o -t- A ))] 

[K4] = 0 

[A^2^] = £'[cosA,cosA.,{(/w7„cos(^o +^,)-mCsin(^G +<^l)) 

-{mT]„ sin(^G +4) + cos(^G +‘f>k)) }((A cos^, - sin6>j + 1)){<I» 

+cosA^cosAjw7„sin(^G +<f>k)-f'<m^^{^G+<Pk)]{^c]^l^ ^ 

[^22] = £'[cosA,cosA„ { mT ]„ sw { ec ., +<^;^)+/w<'„cos( 0 g - t-^) } 

((^,cos0; -/?^sin^, +l)){<Dj{Oj'']£ir 
[^23] = £' )+^^« «>s(^G +<4) }cos A^ {<!> J{^,f 
- sin A, sin A^ {m7„5//7 (^c cosf^f, + 4)}(<1>;}{0 

-cosA,sin“ a4w7„-sw(^g +4)-^-'”Ccos(^o -'-4)}{^'cr}{^9f 
-cosA,cosA^(/..^ -/,^)(sin^ (^c +A)-cos^(^g ->- 


M = o 

[^ 3 i] = £' (cos A, cos A„ [{/^^ cos(^c sin(^^, +<z>,)} 

(- sin (^G + tf )^ ) cos A^ + cos {0^+ (f>k ) sin A, sin A„ ) + 1-/^^ cos{Oa + 4 ) 
sin {Oq + )} {- cos A, cos ) - sin ) sin A, sin A„ )] 

-cos' A, cosA, \{l,^ -/,4cos(^c +A)sin(^c 
(cos' {0^ +<f>,)- sin' {0^ + (ix 

-sinA^ [fe sin (^G + A ) + cos (<9 g + (f )^ )) (sin [0^ + <f )^ ) sin A, 

+ cos A, sin A, cos (^g + 4 )) + {/;;f sin (^g + ) + hr, cos(^g + A )} 

(cos (0Q +<f>,^) sin Aj - cos A^ sin A^ sin (^o + A ))1 ~ (cos' A, cos A„ 

+sin' A, cos A, ) sin' (<9^ + 4 ) + cos' ( 4 , + 4 ) 

+ 2 /^^ sin (^G + 4 ) cos {0Q + 4 )} “ cos A^ cos A„ [{'(f sin (4; + 4) 

+h? +4)}(- cosAjSin(^G +4) + sinAj sinA„ cos(^g +4)) 

+ [hn <^os (^G + 4 ) + hi sin (4j + 4 )} (cos A, cos {0^ + 4 ) 

- sin A, sin A^ sin {0^ + 4 ))] {O, ) {o; f dx 

[MS] = |;a)sA,[{/,, sin (^G + 4 ) + he cos {0G + 4 )} 

(sin ((9 g + 4 ) sin A, + cos(^g + 4 ) cos A, sin A„) + { 4 . sin (^g + 4 ) 

+ 4-7 cos {0G + 4 )} (cos {0Q + 4 ) sin A, - sin ( 4 , + 4 ) cos A, sin A„ )] {<D, } 

[-^ 3 ^ 4 ] = \l \_-m cos A, cos A„ + cos {0a + 4) - sin ( 4 , + 4 )} 

(cos A^ - sin A, sin A^ ) - {mri„ sin {0a + 4 ) “ cos (^g + 4 )} 

cos A^ cos A^ 

["«]=£ /n[cos' A^ cos A„ +sin' A, cos A„ - cos' A^ sin A„]|<I>g}(<l)p}^ dx 

[M 43 J = jj - cos A^ cos 

["“]=£' m cos' A j sin A„ cos A„ + sin A , sin A„ cos A^]{<I)^|(<I>^| dx 



6.2.3 Centrifugal Stiffness Matrix [A'], 


Kn ] = t' [-w(cos- cos^ 0, +COS- 0, cos' A, sin' A., 

- cos 0; sin- A, sin A„ cos j3, cos 0j - 0^ sin 0 ; )) 

+(cos^<9, sin^ A^cos^ A„+cos"^; cos" A^ cos" A^) 

({lm^,cos^(^„+^zi,) + Im^^sm-(^^ + 4)}-{/^_^sin2(^^+«^,)})(0:}{0:}"]A 
[■^12 ] = If [w(cos“ 0J sin A, sin A„ cos A„ - cos0, sin A, cos" A„ 

i^p+^d cos 0j - 0, sin {<E); f + ( COS' 0j sin A„ cos A^ cos A^ 

WmCOs(^G +A)-"<n,Sin(^G + 

-(cos' 0j cos' A, cos' A,)({/^^ -/^^}sin(^G +A)cos{0a+(^,) 
-/,,{cos2(6>„+<25,)}){<t;){(D;f]^ 

Kn] = 0 

[^,4 ] = |!' [w(cos' 0, sin A, cos A, sin' A, ) {<D^ } } 

+(cos' 0^ sin' A, cos A^ -cos' 0, cos' A^ cos' A^) 

+{w7„cos(^o +<z>J-wCsin(^G +<z>,)}{o;}{<l>J^ dx 

[at,,] = I' j^-/w(cos' 0, sin' A„ sin' A^ +cos' 0j sin A^ sin A,){<1>^}{<1>^}^ 

+(cos' 0j sin' A, cos' A„ + cos' 0, cos' A„ cos' A,) 

({im^r sin'(^o +j^^)+Im,^cos'(^o +4>k)] + [Jnc sin2(^o ]fi&r 

[^23] = 0 

] = I"' ( cos' 0j cos A^ cos A^ sin A^ - cos 0j cos A, cos' A^ (>0^ + 0^ cos 0j - 0^ sin 0^ )) 

+(cos' 01 sin' A^ cos A„ - cos' 0j cos' A^ cos' A„ -cos0j cos' A^ sin A^ cos A^ 
{0p+0^ cos 0j - 0, sin 0 , )) {mr]^ sin {0^ + 4 ) + mC„ cos {0^ + (f)^ )} {O^} {<D^ f dx 

[^31] = £' [W„ sin(^G + cos{0o +^,t)}(cos' 0j cos' A, +acos^^ 

+ (e, + ^2 ) { A cos 0j + 0^ sin 0j]-e^ 0^ sin 0j + cos 0j cos A, cos A„ {-a{0, cos 0j + 0^ sin 0j ) 


+e, +e^ +24 +x, y{<I>jKr -[(2cos'^, sinA,cosA,cosA,)({4,-/^4 
s\n{0Q+^0)cos{0a+<^k)+Jr,i {cos^ {0g+^^)-^iv0 {0^ + <4 )}) “ {{^cc cos' (0g + A) 



+ 1 ^ 17 '? (^o + A )} +{2/^c +^^.)sin(^^ +^^;^)|)(cos^ 6j cos^ A^ sin A„ cos A^ 

-cos- 9, cos- A, cosA,(;0^ cos9, - /?, sin^,))]{o ]^& 

[^^ 2 ] = |!' [(cos- sin A, sin A, sin cos(^« + <!>,)] 

+cos~e, sin A,cosA,sinA4/n;7„cos(^o +<^^0-<«sin(^c+(z5,)}){oJ(O,f 

-[(cos^, sin A, cos" A, cos A,)({/^^ -/^^}sin(^^ +A) 

+^7f {cos' {9a + A ) - sin" {9a + A )}) ^-({imcr cos" {9a +A) + 8^““ {^g + A)} 

+ ( 2 / ,jc cos {9a + A ) sin {9a + A )}) (cos" 9, cos A^ sin A^ cos A^ 

-cos- 9, cos" A^ cosA,(/?p + cosOj ->!?,sin0;))]{oj{o;f dx 

[AT,,] = J' [(cos" 9, cos" A^ - cos" 9j sin" A, sin" A^) 

({lm<rc cos" {9a+A) + sin" {9a + A )} " (inicc sin' {9a + A) + cos" {9a + A )}) 

K] = I' [(cos" 9, sin A, cos A, +cos" 9, cos A, cos A„ sin A^) 

{rnr]„ sin {9a+A)+ cos{9a + ^*. )} + (cos" 9; cos" A^ cos A^ sin A„) 

{mrj„ cos {9a +A)- f”Cn, sin {9a + A )}] } {^, dx 
[A: 4 ,] = ||'[w(cos" 9, sin A, cos A, sin" A,){0^}{Oj^ 

+(cos"^; sin" A^cosA„ -cos^ 9, cos" A^cos" A^) 

{w7„ cos(^„ +A)-f^Cr, sin(^r, 

[K^^] = ^'‘^\^m{cos- 9, cosA^cosA^ sinA^){<I>,}{04^ +(cos"^; sin" A,cosA„ 

- cos" 9, cos" A , cos" A^ ) \mrj„ sin {9a +A)+ cos {9a + A )}{‘I ’9 } (^1 f^dx 

K>0 

[A:^^] = [-/w(cos" 9j cos" A, +cos" 9^ sin" A, sin" A,){a>,}{<I> 


6.2.4 Vector : 

L Jl 4 ] 


= Jo' [“ sin cos {e^ + <^, )} COS^ e, cos^ A, cos^ A„ (O J 

m{a(2cos^j -sin^;)cos A„cos^; -(e, +e,) cos^ A^ {fi^cosO; +>!?jsin^^)cos^^^|{0^} 
-{^Vm sin(^G cos(^g +^Jt.)}cosA^ {^c}^{^^m cos((9g + sin(6’G +4)} 




COS A, COS A, 


Z(4),+^. {K} 


\ 7=1 


dx 




-m 


n-\ A 

+ «2 + Z! {h ), + Xk cos^ \{Pp+ Pd cos Oj - /3^ sin <9^ ) (O I 

;=i ; 


|/w£? (2 cos (9; - sin 0 , ) sin (9; cos A^ {<I> J + \mr]^ sin ((9 g + 4 ) 


. ^ "-1 


+/W^„ cos (^G + 4 )} z ih ). + cos Op cos A, {o; } 

V 7=1 


dx 


] = Jo' [{(^ - Irjn ) sin + 4 ) cos (^G + 4 ) cos' 0, cos' A, cos' A, 

{sin' {0G + 4 ) - cos' (^G + 4 )} cos' 0j cos A, cos A, + sin ((9 g + 4 ) 


+/^,cos(^G +4)}sin^; cosA„]{oJfi6f 




W~1 


/^n-1 


-m + ^2 + Z i^e), + + (A cos^, + sin ^/ )|^ Z(^^). + 


COS0; + /?^sin^j)}cos^; cosAjCOsA^jO^j dx 


Using small angle approximation for 4 , the vector ] can be written in a modified form, 
which is given in Appendix A. 



6.2.5 Mass Matrix [a/'] : 

’c 

[m\, ] = I [cos A, {sin Y, cosS, -cos^ (^, cosff, + sin e, - A, sin e,)]- 

sin A. cosAjcosf., + sin f., (^, cos«, + ^ sintf,))-sin A sin A 
{sini^,sin«, +cosn(^,sm», -ftcosff, cose,)}]{®>i 

[M,',] = J[sinA,cosA,{cos(i^,(^,cos«,+^^sin«,)-sin(/,}-cosA,{cos(>',cos6i, 
+ sinf.,(^,cos«,+Asintf,-A,sin«,)}+s,nA,sinAjcosr,s,nff,- 

^in!i^,(^,sin«, cos«, - A, cos«, M*. Jrfsr 
[■^13] = J- [cos sin + / 3 ^ cos sin + sin 0 ^ cos A^ + cos sin A 

t) ' ^ s 1 s 

sinAj{0 

\_^ 2 \ ] = J [sin A^ |cos y/^ + sin y/^ cos 0, + fi^si } t0j )} + cos A^ {cos yr^ ( cos 0, - 

A sin 0, + COS0J ) - sin y/^ sin 0, }j((I)^ ]tAr 

[^22]“/ [sin A^ |sin^^, -cos{/jj. cos^^ +/?^ sin^^ )}+cosA^ {cos^i/’^^. sin^*^ 

-smy/^ {p^ sin 0^ - cos 0, - cos 0, )}]{o ^ ]dx 
L'^23 j = j [sin A^ cos 0, - P^sm0j}- cos 0j cos A^ ] {O J cir 

(J 

[m'] = o 
[wj,]=o 
K]=o 

[M{, ] = J[sin a, cos y^^ {cos y/„ (>9, cos 0j + sin 0i)-pj, sin - sin < 4 ^^ cos } - 

0 

cos A^ cos A Jcos? 4 ^^ + sin y/^ {p^ cos0j + p^ sin6>/ )}+ cos A^ sin A^ 

(cos?i/j,. {p^ cos0f - p^ sin 0^ + p^ cosOj )- sin y/^ sin 0, }](<I»^ )ix 



h 

[A /42 ] - J [sin [cos y/^. cos 6j + sin y/k {Ps <x>sOj + sin 6j ~ p sin OA- cos A cos A 

0 F J 

{sin y/^ - cosw, {j3^ cos 61, + j3, sinO , )}+ cos A^ sin A^ (cos^^. sin 0, 

+ sin yr, {p, COS0: - p^ sinOj + cos$j )}]{<D^ }±c 

[^Is ] = J [sin sin A, - cos A, cos A^ {sin y/^ - cos y/,^ ( p^ cos 0j + p^ sin 9j )} - 

0 S } 

cosOj cosA, sinA^]{<I)^ jtZx: 

6.2.6 Mass Matrix [m^] : 

K]=J[^- cosy/^i^lPAcosOj sin A^ sin A^ -sin^^ cosA^)+ P^{sin6j sinA^ sin A^ 

U 

+ cos cos A^ ) } - 24 sin A^ sin A^ {sin Yk iPa sin 9j P^cosO A - 
sin Yk sin 9^ sin A^ cos A„ (/?^ cos 9^ + P^ sin 9^)+ iPd cos 9j - P^ sin 9A{ 

9, smYk co&9i cos A^, sin A^ -24 cosA^ (Xi%Yk “4 cos A^, sin A^ }]{<l>^}<ie 

[m,', J = J 1^4 sin Yk iPs (cos 4 sin A^ sin A^ - sin 4 cos \) + P^ (sin 4 sin A^ sin A^ 

t) 

+cos 4 cos AJ }-24 sin Aj sin A^ {pnYkiP, cos 4 + y 4 sin 4)+4 sin A, cosA„ 

(>4 cos4 + Pd sin4)cos?i^;(. sin 9j +{pj cos4 -Ps sin4){-4 cosv'*, cos4 cos A^ sin A^ 
-29 j cos A^ sin Yk + 4 cos Yk cos A„ sin A^ )]{«.}<* 

(',) 

K]=JK sin 4 { sin A^ sin A^ + cos A^ sin A^ 04? cos 4 - Ps sin 4 )) + 

0 

(4 ) cos 9j [cos A^ - cos A^ sin A_, 04 cos 4 + >4 sin 4 )] + (4 ) cos A^ sin A^ 

(y^^, sin 4 + y4 cos 4 ) + 4 sin 4 sin A^ sin A^ - 4 cos 4 cos A J{€>^ } 

(O 

[m^,] = J [( 4 )cos?i^;t (>4 cos 4 cosA, sin 4 cosAJ- 2 ( 4 )[cosy^^(; 4 cos 4 + 

0 

pj sin ^ ] cos A, + (sin Yk cos4 )4 sin A,(y^, sin 4 - cos 4 ) + 4 sin A^, 
sin4 [sin A„ (A cos 4 +/?^ sin 4 )]+(4)sin<^^fc(y9^cos4 -^,sin4){<DjtA: 



(O 

[M,^] = J [Oj cos,e, cos +p^ sin^^ cos A J- 2(4) cos cos4 + 

0 

Af sin ^ - 4 cos y/^ sin A, [p^ sin 24 - p^ cos 20 j]+ 

4 sin A^, cos y^,(P^ sin 4 - 4 ^ cos 24 )]iA: 

a) 

[M,\] = J -[4 sin4 cos A, +4 sin A,(4J+4 sin A,(4^ sin4 + >4 cos4) 

0 

4 sin 4 cos A^ 

[^ 32 ] = 0 

K]=o 

[m;,]=o 

[^M 4 , j = J [4 ^os Yk [ (cos 4 cos A^ sin A^ + sin 4 sin A^ ) + 4# (sin 4 cos A^ sin A^ 

0 

- cos 4 cos A^ )] - 24 cos A^ sin A^ cos Yk iPd sin 4 + A cos 4 ) + 

4 sin cos 4 cos A^ cos A^ (4^ cos 4 - A sin 4 ) ~ ^4 cos sin A^ (4* sin 4 

-4rf cos 4 ) - 4 (sin Yk sin 4 ) cos A^ cos A^ (4, cos 4 + Pd sin 4 ) + 
cos A J cos A„4 sin y/^ (4j sin 4 - Pd cos 4 )] {*1*, } ^ 

(O 

K]=J[^, sin Yk [ P s (cos 4 cos A^ sin A^ + sin 4 sin AJ + 4^ (sin 4 cos A^ sin A^ 

i) 

- cos 4 sin A^ )] + 4 cos A^ cos A^ cos cos 4 (At sin 4 “ Pd cos 4 ) “ 

24 cos A^ sin sin A^ (4s cos 4 + Pd sin 4 ) "^ 4 cos Yk sin 4 cos A^ cos A^ 

(4, cos4 + 4^ sin4)- 24 sinA, svciYkiPs sin4 “Ad cos4) + 4 cos,Yk cos4 
cosA^ cos A, (Ad cos 4 “As sin 4 )} 

r ' W\ 

[a /^3 ] = j [-(4 ) sin 4 [cos A, cos A^ (Ad cos 4 - A. sin 4 ) + cos A, sin A J - (4 J cos 4 

0 

[sin A , + cos A J cos A^ (As cos 4 + Pd sin 4 )] + cos A^ cos A 34 ( A^ cos 4 + 

AdSin4) + 4 cos4 sinAj+4sin4cosAjSinAa]][Og}<3&: 



6.2.7 Mass Matrix 



i, 

\j^n ] = J {( [sin Wk sin A, sin (sin 9, sin A^ - cos^, cos A^ cos A J - 

0 

cos A^ (cos A^ cos^jj. sin Ag -sin^^ sin^^^ cosA^ cosA^) + sin A^ sin A^ 
(sin y/, cos 6, cos A, cos A, - cos y/, sin A, )] + (a) [sin 6^ cos y/^ cos A^ + 
cos 9j cos y/^. sin A^ sin A ^ ] + (e^ + (w), ) [sin y/^ sin A^ cos A^ + 

sin y/^ sin A^ sin A^ cos 9 j ]} (<1>^ }£& 

K 

[M ,^2 ] = J {( xJ [cos cos A^ sin A^ (cos 9^ cos A^ sin A^ - sin sin A J - 

0 


cos" A^ (cos A^ cos y/^ sin 9^ + sin A^ sin - sin A^ sin A^ 

(sin^j, sin A^ +cosy/^ cos 9j cosA^ cosA^)] + (a)[sini9j siny/'i^ cosA^ + 
cos 9j sin sin A^ sin A ^ ] - (e, +e^+^ (u\ ) [cos y/^ sin 9, cos A^ + 
cos y/,^. cos 9j sin A^ sin A ^ ]| (<I>^ } dx 
h 

[a/ji ] = J (( ) [sin sin A^ (cos 9, cos A^ sin A^ - sin sin A J - 

0 

cos A^ (sin y/'^. cos 9j cos A^ cos A^ - cos y/,^ sin A^ )] + (a) cos 9j cos y/,^ cos A^ 
+(e, + ^2 + 2 (“). ) sin Wk ^os cos A ^ 

] = J {( ) [cos y/,^ sin A^ (sin 9j sin A^ - cos 9j cos A^ sin A^) - 

0 

cos A^ (sin y/^ sin A^ - cos y/^ cos 9j cos A^ cos A^ )] + (a) cos 9j cos A^ sin yr^. 
-(e, + ^2 + 2 (")< ) 6>f cos A J {O J iA 

[a/^j J = J|-( ) [sin A^ (sin 9j cos A^ sin A^ + cos 9j sin A^) + cos A^[sin + 

0 

cos sin ^ , ) + cos A,cosA3sin^j]]+(a)sinAg- 

(c, +Cj +5](K),)sin^; cosAjf®,}* 

[«3.] = 0 

[Mi] = 0 
[Mi] = 0 



[■^41 J ~ J{( ^4 ) Vk Aj cos (sin 0j sin A^ — cos 9j cos A^ sin A^) + 

II 


cos A, sin A^ (cos y/^ sin A^ - sin 0^ sin y/^ cos A J + cos A, sin A„ 

(sin y/, cos 0, cos A, cos A^, - cos y/, sin A, )] - (a) [sin 0j cos y/, sin A, - 
cos 0, cos y/^ sin A„ cos A J - (e, + e, + (w), ) [sin y/^ sin 0j sin A, - 

sin cos A j sin A^ cos ^ ; ]| 

h 

[Ml. ] = J {( ^ J [cos cos A^ cos A^ (cos cos \ sin A^ - sin <9^ sin A J + 

0 

cos A^ sin A^ (sin y/^ sin A„ + sin 0^ cos y/^ cos A^ ) - cos A^ sin A^ 

(sin y/y sin A^ + cos y/^ cos cos A^ cos A ^ )] - (a) [sin 0, sin y/^ sin A^ - 
cos 0j sin y/f^. sin A^ cos A J + (e, + e, + ^ (z/), ) [cos sin sin A^ - 
cos Yk cos A^ sin A^ cos ^ ^ | tA 

] = J {( ^^- ) [cos A^ cos A^ (sin 0 ^ cos A^ sin A^ + cos 0 , sin A^ ) + 

0 

cos A^ sin A^ j^sin A^ (/?^ + cos sin ^ - cos cos A^ ] - 

cos A^ sin A„ j^sm ^ + cos sin ^ ; ) + cos cos A„ sin 9i ] - (a) [cos cos A„] - 

{e^+e.+ Zc-) ,){cos0j sin A ^+cosAjSinA^ sin^;]||C)^|(ir 


6.2.8 Matrix [M^] 14x3 i 


If f 

[M,"*, ] = j [-m cos A^ cos A^ {a cos - (>9, cos sin 0,+\) ^ (/^ ), + 

e\ V J~1 


-(e, + 1 ? 2 ) cos Yk + ^ cos A^ cos A^ (/, ), + ^ 




{*.)* 


[a/, j ] = |[-mcosA, cos A^acosO, ~(P,cos6, - fijSm&i + 

0 


1 ) 


fs(o,+^. 


v.=i 


j 


-(e, +e2)sin^i/jt 


[ 

/n-l 


]■ 

/wcos Aj cosA^-^ 


cos Yt ( 

1 

V.=1 

J 

J. 


{04* 




K]= 

[K]- 

[K,} 

[K] 

[<]= 


V 

J [wa cos cos Oj (2 cos (9^ - sin <9^ ) }<ir 

0 

f T""' ^ 

J [-mcos A, cosA„ cos {2a sin cos^^ + cos^i^jfOj^ 

0 V -=i ' y 

f r ^ 

J[msinA^cosA^cosA,(2asin^i^jCOs^;- '^{h),+Xk sin^iAj,.}{<Dj<ir 
0 V 1=1 " y 

/ n~l \ 

j[-wcosA,cos^A, 2(^0, (Pp+J3,cosej-P^sin0j0j)}{^,}dx 


V >=1 


0 

0 

0 

le 


^ n-~\ 

f [-W cos Aj cos A^ cos 0j (a cos cos 9j - sin $^+1) ^ ), + 

0 V!=l 

-{e,+e^)&mil/^[^^]'\dx 

k fn-\ 

f [m cos A^ cos 9j (a cos 9j - (>9^ cos dj - sin 9j+l) ^ (/* 

0 Vi=i 

-(e^+e^)cosy/;,{^^}~\dx 


+ Xt, 


\1=1 


6.2.9 Vector 

/2-I 

r^n ~ f'[-/wcosA^ cosA3[(3a-(e, +C 2 )-((^^ cos^; + sin ) + IX^ (/^ X +^*:) 

L J Jo 1=1 

-{Pj (cos (9; + sin 9j) + Ps (sin 0j - cos ffj ))} cos A,0j 

-{((P, cosO, +P, sm9j)-l)(^(.lel +x,M cos" A, cosA, 

1=1 

+ {a + (P^ COS0J + Pa sm0i)(^(le)i cosA^ cosA^]{<I>^}£A: 



[Vl,] - |,''[-'»[{« + (AcOsS, + (/,X 

/ I 

-{a cos e, cos A, cos A, - {{p^ cos 6, + p^ sin di) + \)[^ (/, ), + x/ 

. V f I y 

“ (^1 + ^2 )}<?/ cos A^ sin A^ sin A^ 

-(Acos^, + A^in«,)(i;(/,), +x.)4=cosA,1{4..}* 

/-I 

[^'] =° 

Jo I- +-*^<:)(A cos^^ sin^^)}(9^ cosA^ cosA^ +aOj 

n-\ 

-{(e, +e3 +Z(/J, +x,)(l+(Acos^, +y^^sin^,))-(e, -^e, -a)dj cos' A, 
+ a sin e, cos 9, 0; cos A cos A J (O,^ }dx 



The integrals of the shape functions, which are needed for the evaluation of the mass 
matrix, are given below; 


dx = 




dx = 

0 

0 




r I ^ 


V60y 


3 


9 


_A. 

35 

60 

140 

70 



A. 


/; 

60 

280 

60 


280 

9 

1. 

2 


_4_ 

140 

60 

7 


28 


/; 

4 


1_ 

70 

280 

28 

168 

" 9 

52 

2 




12/, 

4 



J -9 

88 

68 


Je 

-16/, 

- 64 _ 


' 0 

-24 

-6 ■ 



4 

-4/e 

-2/e 



0 

24 

6 



.-4 

0 

64. 






60 J 


1 

0 

-1 


0 

16 

4 


-1 

4 

7 



■ 156 

22/, 

54 

-13/, 

r 4 1 

22/, 

4/e^ 

13/, 

-3/e^ 

U20J 

54 

13/, 

156 

-22/, 


-13/e 

-3/; 

-22/, 

4/; 


64 






dx 


0 


l3o; 


-30 

-6/, 

-30 

6/. 

0 

-6/. 

30 

6/e 

30 

.-6/. 


6/e 

'4 2 

-1' 


2 16 

2 


-1 2 

4 






11 

20 

-1 

'/el 

/e 

4/e 

0 

60j 

-1 

20 

11 


0 

-^/e 

-h 



f— 1 

UoJ 



-48 

-6 

-4/e 

-3/e 

48 

6 

-4/e 

7/e 


^ 6 
6 
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Figure 6.1 Finite element model of a blade 



Figure 6.2 Element nodal degrees of freedom 
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Chapter 7 


Formulation of Element Matrices Associated 
with Strain Energy Variation 

The elemental matrix associated with the strain energy variation is derived by 
substituting the assumed expressions for the displacement function in the strain energy 
variation SU, and carrying out the integration over the length of the element. The 
resulting variation of the strain energy can be written in the form: 

w. ={<??)' (7.1) 
Where, is the elemental stifl&iess matrix, is the nonlinear stiffiiess 

vector and represents the vector of unknown degree of freedom and is of the form: 

The nodal degrees of freedom are shown in Fig. 6 2 

7.1 Linear Stiffness Matrix 

[^,^] = j[(£^cos^o-^sin^<,){(I):}{0:f]jx 

[ii:f,]=[ii:f,]=j[(^cos^o-:^sin^^){<i>:}{o:}"^ 

0 

0 ^ 

[a: 5] =[A:f,] = j[r(-5^){4>:){4>;r V 

0 





[at' ] £ 1D'){*;J{4,. J' J,_, 

[^S]=[^.‘]=/[(-Sc){®;){4,;f jfc 

K}{4.;f {4.;}{<i,;}' 

(r=£4£);-GJ){o;){4.;jrJ* 

[^s]=/[(^){<i>;}Kr}fe 


If the properties like El, EI,„. Eljt etc , vao- linearly along the elenrent length; 
quantities defined in the linear stiflfhess matrices can be expressed as; 


El, 








EI„„ = 

m 


+ 


+ 


+ 


El . + 








/ 


COS^Q - 


COS^Q - 




EE^-EI, ^ 

f'h Cn- 


EL,, + 

V^h 




sin^Q + 


sin^Q + 




^ n.-EE ^ 

‘^’h ‘,ni 


El 

nti 


EL„ -El 


sin ft. 


sin ft. 


cos ft. 


cos ft. 


EAD^ = 




EAD[. 


EAjj^ = 


EAD{,+ 


( EADl,-EAD[,^ 


cos ft. 


cos ft. 


£4 A, +^E^2i-EAD^^'\ 


EAD'^,+ 


( EAD'^-EAD', 
I 


sin ft. 


sin ft. 


+f ~" . '^•' 


'I 


COS ft. 


EAt}^^ + 


EA7 ]^^-EAtj„, 

I 


sin ft. 


( 7 . 2 ) 

various 




■ 0 2 0 - 2 ' 
f \ \ -2 -\ 2 -1 

0 -2 0 2 
2 1 -2 1 _ 


CC\ 





0 



12 

6 /. 

-12 

6 /. 

rn 

6 /. 

4 /; 

- 6 /. 

2 /; 

22 

-12 

- 6 /. 

12 

- 6 /. 


.2 

2 /; 

- 6 /. 

4 /; 




60 


36 

6 /. 

-36 

0 

6 /. 

2 /; 

- 6 /. 


-36 

- 6 /, 

36 

0 

0 


0 






-18 

6 /. 

18 

- 24 /,' 

1 ^ 

9 /. 

- 2 /; 

9 /. 

-7/; 

130/ J 

18 

- 6 /, 

-18 

24 /, 


. 21 /, 

8 /; 

- 21 /, 

13 /,' _ 


]x{^:mfdx: 




6 

2/, 

-6 

/e 

2/, 

/; 

-2/e 

/e 


-2/, 

6 

-^/e 

.4/. 


-^/e 

3/e^ 


dx= 


^ 2 ^ 


3 

-/ 


0 -3 

2 /. -K 

0 3 

-2L i 



0 


22 


-4 

0 

4 

0 


8 

0 

-8 

0 


-4 

0 

4 

0 


70 







4 -8 4 

3 ^ - 4 /. /. 

-4 8 -4 

K - 4 /„ 3 / 


l{^:}{o;fdx= 




0 0 0 

- 4 /. K -41, 


0 0 


0 


4 i - 8 / 4 / 


|*{®:)K}*= 


J_ 

60 


18 


9i 




/ 1 ^ 


v.3^, / 


24 

- 42 ' 

12 /. 

- 11 /, 

-24 

42 

- 28 /, 

19 /. J 

12 

- 6 ‘ 

2 /. 

-/. 

-12 

6 

- 2 /. 

/. 


H®:}{®;}’'*= 


^ 1 




3 -12 9 

2 /. - 4 /, 2 /, 

-3 12 -9 

/. - 8 / 7 / 






4 

0 

-4 

4 /, 


-8 4 

0 0 
8 -4 

■ 8 A 4 / 
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IKKK 


0 


_L 

3 / 


7 

-8 

1 


-8 1 ■ 
16 -8 
-8 7 





8 

0 

-8 


-4 

0 

4 


}{<!.; }{o;r,y.=l 

0 

0 ^ 

0 


3 

-4 

1 

-4 

16 

-12 

1 

-12 

11 


16 

-32 

16 ■ 

-32 

64 

-32 

16 

-32 

16 




fl ] 

■ 8 

-16 

8 ■ 

-16 

32 

-16 

UJ 

i 

8 

-16 

8 




-2 -8 

10 

4 16 

-20 

-2 -8 

10 
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During integration along the length the cross-sectional const 

constants were assumed to 

vary linearly along the element Using the linear part of the en . i . r xr. ^ . 

, . - , 1 . . term), and 

replacing the vanables with their corresponding nodal vain. , . • , 

, „ . aiues multiplied with the 

appropriate shape functions the entries of the elemental stiffn.. 

ness matnx were found. The 

matrix entries for J were found out to be the following' 


= [eI^^ cos0a -EI^, sin^^)[^Hj+ cos^<, -^sin6»3)^4j 

= ^3.1 =K,cos^^--:^s.n^^)[^-^j+A(^c^^ 

= =[eI^^ cos^g -£7,^sin^o)[^|-j+ a(^cos^o '^sin^'^jjj^y 

= ^54 =(^,'',cos^G -^77sm^G)(^^j+A(£/^,cos^o-£^sin6'G)^^ 

= -^6x1 ={EIi^coseQ -£/,,sin^G)[^|-j+A(^cos^o -^sin^'o)!^- 
= ^7x, =(^cx)s^c -^sin^«)j^~j + A(aJcos^Q 

= ^s.i =[^ir, cos^G -EIr,r,^mea)[^ + ^[EI^^cose^ '^sin6'G)[^y^ 
f a\ /i'\ _r A\ 

^ +AtoEAD; - +EADj( 0 ) + AEADj ± 

vJ UJ {ly 

r o^ r a\ r X 


" -^10x1 "■ *^0 


— ^ 
= — - Tq £/iDi y 

V I 


,EAD[\ -4VAroM);f-yl+£4D,(0) + A£inr-i. 
t) IJ /2 

'N , m I A ^ f\\ t A^yrR ■( 4 ^ 


f ("4-1 + Ar„£4£l|’( 7 I + £/<£>, (0)+ A£1 d /4 

V J Vv \ 1 ‘ 


= ■^12.1 =-^7a 


'4>| 
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= [eL^ cos^g -EI,. sm^o)(^yj + a(£/,^ cose^- EI^. sm0^){l) 

--^3.2 =[EJ^.co%Oa -■^^<sin^o)(^-|-j+A(£^cos^G-^sm^G)( 
'^4-2- COS <?Q - El sin j j + A ^£7^^ cos - El sin ^^ ) ( 0 

= ^5 - 2 = in cos Oo - El sin j + A (£7,^ cos (% - El sin ^a){j 

= ^6.2=(^cos<9o-;^sin^o)[^yj+A(^cos6»<,-£^sin^c)0) 

= = (£7^,^ cos^G -£7^,sin6'G)|^-|-j + a(77^,^ cos^g -EI,^ sin^G)[ 

■■K^.2= {EEr, cos Oq - El^^ sin 0^ ) fyl + a(£7^^ cos 6^ - £7^^ sin (9 g ) ( l) 


As'i 72^ ( 4^ 

■ K^^_^t^EAD[ - +^r^EAD[ - +EAD^ +A£4£),(0) 
\l ) v3j \ I" ) 


■- = t^EAD'A --UAt^EAD'A -- +£47), 4- + A£4Z), (O) 

\ U V 3j \l- ) 

: r,,., = r„£55;flV 4r„^f|V^f-4'l+ A£4^ (0) 
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-^3x3 = cos 9q - sm 6 ^ ) j + k(eI^ cos 6^ - sin 9 ^ ) j 

■^3x4 ~ -^4x3 “ cos^Q —EI^^ sin^Q j+ cos^^ 

■^3x5 “-^5x3 j + A^£Z^^cos^Q-£/'^^sin^G^|^— 

^ 3 x 6 = ^ 6 x 3 = [eI^„ cos 9 a -EI^ sin^G)J^-|-j+ a(^cos^g 

■^ 3 x 7 = -^ 7 x 3 = cos 9c - EI^^ sin 9 ^ ) j + A {eI^^ cos 9^ - EI^^ sm 9^ ) j 

-^3x8 =-^8x3 + 

( 4 A f I'N f a\ 

^3x9='^9x3=^o^^i -TT +ATo£4D;— - +EAD,{0) + AEADi -T 

V W \ U W j 

-^3x10 ~ -^10x3 “ Tf^ElAD-^ ( ■^1+ AtqEAD^ — + EAD^ (O) + AE4D, — 

\* J \U V* J 

■^3x11 ~ -^11x3 ~ ^0 ^ + Az'q£( 4Z), - + £(4Z)j (O) + A£4Z)[ — -pp 

V ' y \ U V ‘ y 

f 4A f 

■^ 3 x 12 ~ ^ 12 x 3 ~ — AEAT]^ — - 

V ‘ y V 'y 

fS^ ("4^ 

^ 3 x 13 ~ -^ 13 x 3 ~ TT ~7 

V‘ y V'y 

=-^f-4VA^r-y' 

V » y V. »y 
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K, 


■^4x6 
^4x1 
^4x$ ~ 


= [eI,, cos0^ -W;, Sin^^)[^lj+ a(;^cos^^ -:^sin5»^){3) 

= -^5.4 = cos^g -El^r, sin^o)[^|-j+ a(£^cos^g -^sin^'o )j^y 

= i:,.4=(£f^,cos^o-:^sin^^)(^|j+A(£^cos6»,-:^sin6'^)^^^ 

--^7x4 -{eI^jiC^sOq -E/^^sin^g j+A^£7^^cos^c 
= ^8x4 =(e/^,cos0g -^■^7>7sin^G)fjl+ a(1^cos^g -£^sin6»G)(3) 


-1 

’/J 


■^4x9 ~ - 


/'A f a\ /■ . 

K^^=r^EAD{ - +AtoE4D; - +EADA^ +EEAD, - 
VA V.3j \l- ) \I 


■^4x10 “ -^10x4 ~ T + 4TqE1ADI 


\ 


I 


f--l 

3, 


+EAD, 


( ^ 

l- 


V i J 


K, 




4x11 


: = roE4Z),|yJ + AroE4Z); 

-^n .-=rT-f\^ 


^1} 

U; 


+EAD, 


^4^ 

\l'-J 


+ AE4Z), 


+ ^EAD^ 


(- 

[h 


f--l 

I h 


^4x12 ~ -^12x4 


•^4x13 “ 


K. 


4x14 


^13x4 =-EA7J, 

K. 


\h 

(-- 
K h 

^ 3 ^ 


-AEAjj^ 


-AEAij^ 


UJ 

-r_8> 

I 3; 


14X4 = -EAjj, l^yj - AEAjj^ - 
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^5.5 = cos sin 6*^ ) I^H j + a(£/^^ cos^3 + 

■^5x6 =-^6x5 +^f,sin^o)|^— j+A(£7^^cos(9G +-E^^,sin6'o jfy 


^5x7 ^7x5 ( COS 6q + £7^^ sin j — j + a ^£7^^ cos + El sin ^ ^ 

^5<s - ^$-5 - {^^r,n cos Og + EI^^ sin j j + A ^£7^^ cos 0g + EI^^ sin j f yl 


^5a 9 ~“ ^9x5 "“ TqEAD^ ^TqEAD^ l^y + EAD^ (o) + AEAD 2 J 


7( 8 


.f A 

■^5x10 ~-^10x5 ~ | + 7 +-£ 47 ) 2 ( 0 ) + A£4Z)2 ^ 

\ ‘ J \ U \ E 


^5x11 =£iK5 = ro^^2|^yrJ + ^^o£4£2l yJ + £47)2(0) + A£^2h 


£ 5 x 12 = £ 12 x 5 =-£4C 7T -A£4C 


^5x.3=£,3x5=-£4C ~ -A£4C 

\ I J 
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k,.,=[ei, 





cos(9g + £f^^sin6>G)(l) 
^6«7 =^7.6 =K7Cos^o+^sin^^)[^-|-j+A(^cos6>^+:^^ 
^8^6 - {eI^^ cos Oq + sin j j + cos 9^ 


Ke.s=K,. 




^G+-£^f,sin(9o)(l) 

-fS^ f 2 ') 

: - +Ar,EAD'A± 

\U 'UJ 


4 1 f 4> /o 

-WAt^EAD^ — +EAD^ 4 

U V 3y \f 


-( 1^ f 9'\ 

.,„6=roM)' - +AToEAD; - 

\l J y3j 


-EA^a 


f3^ 

u 


J 


-A^.[f) 


I 


J 


m 


= /^.3.6=-£:4c[^-y]-A£4c(^~ 


:K,,^,=-eacJI]-aeac 

V' J 


'(-] 

UJ 


El, 


,,cos6>y 4-£:/^^sm^o)te 


jy^{EI,,cos 9 ^+EI^ 


^nn 


,,EAD'. 


cos 9^ + EI^^ sin ^0 )f -y-j + a(£/, 
/ 


V 


/= 


-EAC, 




+ EAD^{0) + AEAD^ 

eadJo)+aeadS 

\l) ^ \E 

{ 2 ^ 

+ Aro£4Z)2 — +EAD2{0)+AEAD^ 

V u 

n 



,,,jC 0 s 9 a+EI^ 


J 


J 


E4cf-|l-A£<. 

V ' 


J 

( 4 \ 

~ -aeac„ 

V I J 


I /, 

'(-] 
UJ 
/ ^ 


I /. 
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w 4 ^ 

EI,, cosOa +£/^,sin^o)U U a(£7,^cos^g +^sin6'o)(3) 


■^8x9 “ -^9x8 


t^EAD’^ 


rn 

UJ 


+ h.TQEAD[ 


•^8x10 ■“ -^10x8 “ '^oEAD2 


I 


■\-l^tnEADl, 


rn — 

- +EAD, 

UJ 

^ 8^ 




+ LEAD. 


K.^u =^iK« =t^EAD' 


E-SxU ~ -^12x8 




f3^ 

U 


+ AtqEAD!, 


V 3 

U. 


+ EAD.^ 


r a\ 


I"* 

■/- 


+ AE4D-, 


+ EAD’, 


\Ej 


+ AEAD. 


i^8xl3=^.3x8=-£^C 


v/y 

^ 4^ 


A£4C 


n 

39 


1) 


a:3.„=^hx8=-^C 


V 

f-1 

l/J 


fl\ 


■aeac. 




U] 

I u 



K,^,={GJ + ^EAD',)\1 \ + a{GJ + E^EAD;)\ E \+t,EAD, 

Z._Al 

3/ 3/J 


i._d 

I- P 


+ 


f 2 2 

Ar^EADJ — +EADJ ^ | + AE4Z)J ^ 




+ TqEAD^ 




AVqEAD^ 


u 

Gil 3 / ) 




^3^ + A£4D3| 


16 

P 



4 4 

+ r^EAD^\ -p+JI 


+ 


Z' n\ 


Egx\2 ~ -^12x9 ~ “^qEIADq 


^9X13 =^.3X9 =-^0^0 


•^9x14 ~Ar,4^p — TqEADq 


v3/y 

3/9 


-4t„£ 10;| i 


-£:4Z)„ 




-A£4Z)n 


-Aro£^; -- -EADo{ 0 )-AEAD, 
V 3 y 


I 3/y 
3/ 


2 

3/ 


'N 

rn 


r4^ 


fio^ 

-ATqE4Dq 


\-EAD^ 


j-A£4Z)o 


J 

^69 


U"; 


l3/ ) 
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K, 


10x10 ^10-10 i 


ro £4^; + A(a/ + r,-£4D3' )f|V t,,EAD, (O + O) 


+ATcjEAD^ 


16 16^ ^64^ f^~>\ 





16 

’/- 


■^10x12 “ -^12.10 ~ -t^EAD'^ 


3 h 


■Ar,EAD', 


(--] 

3J 


-£4Z)„ 


v/V 


-AEAD, 


r a\ 


^IJ 


^mu ~ -^13x10 ~ '^qEADq ( — 

V3/y 


-AtqEADq 


/'c^ 




-EAD^{0)-AEAD,, 


f-1 

1 3/. 


•^10x14 ” ■^i4'io “ ^qEADq 


f Q ^ 


\ 31 j 


■ ^t^EAD[{-l)-EAD,, 


f fi ^ 


V ' ) 


-MADr, 


’3/ 


1x1 1 ■“ -^1 ivi 1 {GJ + TqEAD!^ j — 4- a {GJ + TqEAD^ + TqEAD^ 


.3/ 



4 4 

r-^r- 


31 31) 




■^11x12 “ -^nKii 


~ro£4Do 


fl] 

-Aro£4£>o 

fl] 

-EAD^\ 

f 4^ 

-AEADJ-— 

l3/j 




V EJ 

1 31) 


^11x13 "“^13x11 “ -TqEADo 


£1 ivi d — Kj / 


■-T.EAD' 


r 

X 3/y 
r j\ 


- At^EAD;, (-2) - EAD^ (0) - AEAD^ 


^3/y 


-AtoEAD' 


411 


■EAD, 


r a\ 


\i-j 


-AEAD^ 


31, 


\3Ij 


K, 


12x12 


■EA 


V 

\3l j 


+ AEA 


r^\ 




■^12x13 ~ -^13x12 ~ EA 


. 3/, 


r 


+ AEA 


3) 
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-^12x14 ~ -^14x12 ~ 


r 1 \ 


v3/; 


+ AE4 


.6. 


-^13x13 “ 


^16\ /8^ 

+ A£4 


v3/y 






13x14 '‘^14x13 


K,.,,,=EA 


31 


+ AEA{-2) 


■^14x14 ~ 




31J 




(7.5) 


7.2 Nonlinear Stiffness Vector 

The nonlinear stifihess vector is given by the following sub-vectors: 



{^'} = J j{v. + 0.5(v; +i-;)) -4E4D, + 0.5,g£4C„- 

->.'.(£44:+ - v„ 

+ sin(9(3 (v„ cos^G +^» sin^G)[<^x -GJi)-(I>qGJq'^{<1>',] 


+ |cos 0Q cos 0Q - sin do)[UGJo+GJ,)+<l>,GJ,'\ 

-0.5(f^ + m>1')EAtj^ -0.5^^^F4C, +^[-<^„F4£)2 + («:, + 0.5(v^ +wl'fjEA^^ 

+ 0.5il>lE^,-w„ [(^+^)sin^G +(^-^)cos^g] 
-v„[(£?;;-£^)sin^G+(^+^)cos^G]}K}}^^ 
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} =}{{'''• + + 0.54>lEAC, -4,,z,BADl 

-w^ [eA^^ + (pEArj^ ) - - <fEA^^ 

+ cos Oq (v^ cos Oq + sin 0 ^ ) \j>^ (GJ^ +GJ,)+ <j>^GJ^ ]} } 

+ {sin^G K cos^G -V, sin^G)[^^ (GJo +GJ,)^(j>pj^'\ 

-0 5(v^" + w;)£4C + + 0.5(v^ 

-<PJc>EAD{ + O.S<l>l'EiC,-wJ^^^+W^^ +(-£J,^+^)sm^G] 

-^J\[^Ics -(£^+^)sin^G|K}}^^ 


le 

{^3^} = l{v«[-^„^+(«, + 0.5(v^+w,^))liC-^5,ro£:i^+0.5^5,^:^] 
0 ^ 

-«'xx sin^G +£■/,, cos^g]-v„[£^^^ sin^c +EI^^ cos^g] 

+(l>[-<l>^EAD,^[u,+0.5{vl+wl^^^ 


-wJ{EI^^ +£7^^)cos^g +i-EI„r, +EI^^)sm9a'\ 

-EI^^)cos0a -(EI^^+EI^^)sm0^] 

+("x +0-5(v^ +^1))EAti^ -( f)^r^EAD[ + 0.5(f)lEAC^ 

+«'« cos ^G - El nr, sin \eI^^ cos 6>g - El^^ sin 0^ ] 

+^\^4>xxEA.D2 +(«x +0.5(vj +M'^))£44'a -(f>tEAD2 + 0.5(f>lEAC2^ 

->*'« [(^ + ^) cos 0o + [-EI^^ + )sin ^G ] 

-Vxr \[eI^^ -17^)cos^g -{EI?n +^-;f )sin^G]}{^,} 

+ |-r„ j^O.SE^o (^x + ^x ) + 0.5EAD\^l - <l>{w^EAD[-v„EAD’^ 

+<t>^ \^-(f>^EAD^ + (i/, + 0.5 (v^ + wl ))£4Co 

-4>xXoEAD\ + 0.5(I>IEAC2 -w„ (eAC^-^ ( j>EAC^-v^ [ eAC, - (j>EAC2)^{<l>n) 
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(7.7) 


+ {-0 {k + k )-0 -v„£tD,)}{^5;}tfe 

{^/} = I {0-5£^(v^' +w:) + 0 5 (P^EACq 

The underlined terms in the above sub-vectors are those, which are associated 
with the axial strain (w;, +0.5(v^^ + vi';))at the elastic axis of the blade. These nonlinear 

terms are modified to a set of linear terms by substituting the axial stress by the axial 
inertial force. The description of the substitution approach is given in the next section. 


7.3 Linearization of Nonlinear Terms Associated With The Axial Strain 
At The Elastic Axis 


The equation of motion corresponding to the axial degree of freedom can be wntten 
symbolically as; 


^0 

mQ^ 


[-kI-z.-o 


(7.8) 


The axial stress resultant can be wntten as 


i;=EA[K,+0.5(v’+w=)]+f (7.9) 

Where F contains all the additional higher order terms which can be neglected. The 
simplified expression for can be written as 

jF « Ea[w^ + 0.5(v^ + w; )] (^ lO) 

Neglecting all the higher order terms and also the time derivative terms, the distributed 
inertia force can be written as: 
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where 


^ 1 “ ~ Pd sin 0 , )} cos 0,' cos A^ cos A„ 

Q ~{(^i +^2)“<^(A +y?jSin^j)^cos^,"cosAjCosA^ 


Where the symbols have the usual meanings.(are given in chapter 5) 

Integrating eq.7.8. 

i 

F= \Z^dx (7.12) 

From the above we see that can be written as; 

V=f{^) = a,i'-+a,^ + a, (7.13) 


where. 



a,, <3,, a, are given as; 


mQ.' 


a=- 


E.. 


2 


mQ' 


^2 


-cXiO -Q 


l=l 


(7.14) 


a. 


mCZ 



[i (o,- 1 S(0] 



so finally we get the axial strain terms as a quadratic polynomial: 


k+0.5(v'+w^)]=^ 

EA 


( 7 . 15 ) 


Now linearize the underlined terms using the above relation to substitute the axial 
strain term. The following additional integrals apart from those already listed earlier(6.8) 
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35 / 
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-1 

14 

-T 
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-/ 
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70 
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1 
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1 
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■ 6 

1 

-6 

1 

5/ 

10 

5/ 

10 

1 

2/ 

-1 

-/ 

10 

15 

10 

30 

-6 

-1 

6 

-1 

5/ 

10 

5/ 

10 

1 

-/ 

-1 

2/ 

.10 

30 

10 

15 

■ 1 

—7 

> 

11 


15/ 5/ 15/ 

-2 32 -26 

15/ 15/ IsT 

j_ 

5/ 15/ 



0 


1 

-2 

1 

2/ 

3/ 

6/ 

-2 

8 

-2 

3/ 

3/ 

/ 

1 

-2 

11 

6/ 

/ 

6/ 




7 

-8 

1 

3/ 

3/ 

3/ 

-8 

16 

-8 

3/ 

3/ 

3/ 

1 

-8 

7 

3/ 

3/ 

3/ 


(7.16) 


From the underlined portion of the nonlinear stiffaess vector, following was obtained: 
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I 

From the above we get the following additional terms which to be added to the elemental 
stifi&iess matrix; 


^1.1 = 


12 

35/ 


3 

5/ 


•a, H a, H CL 

I Cl 2. jmj 3 


6 

SI 


K, 


U2 


K. 


1 , 1 1 
■K^^y=—a,+ — a + — a, 

14 ‘ 10 ' 10 ' 

P. 12 3 6 

— o.^ 


1x3 -"’"3x1 


35f * 5/ 


SI 


Ku 4=K. 


1 


1 


4x1 - - 


K^.2 = 


21 


K. 


2^3 


35 ' 10 ^ 

h 21, 

— ^a, + —a, + — 

105/^ 30 15 

1 1 1 

14 * 10 " 10 ' 


K. 


I 


I 


2x4 


1 

30 


^3x3 


■ '^4x7 — ^1 ^7 

70 ’ 60 ^ ^ 

12 3 6 

-a, + — a, H a. 


357 ' 57 


SL 


35 1 3 


‘‘4x4 


jr -34 4 2/, 

Ci-f H d'. 


35 ' 10 ^ 15 ' 
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■^5x5 "^1x1 


6x5 ■^1x2 


K, 


5x7 


^ 7 x 5 =^- 


1x3 




8x5 — ’ -^1x4 


Ks. 6 =K-. 


K^,=K. 


2x2 

7x6 “ -^2x3 


ii: 6 x 8 =^l 


8x6 -^2x4 






“■Sx? 


■^ 8 x 8 "^ 4 x 4 


^ 9 x 9 “ 




£4 


3x4 


1 1 7 

~< 3 (, H 61 , ^ 61 , 

5 /, 3 /, 3 /, ^ 


_ _EAC, 

■^ 9 x 10 ~ -^KVO — 


*-10x9 


•^ 9 x 11 ~-^nx 9 ” 


EA 

EAC. 




57 


37 


-a-. 


3 /. 

1 

, -a, + — a, + — a, 

5 /, 6 /, 3 /, ' 


-a. 


1 


1 


■^ 10 x 10 ~ 


EAC, 


EA 


EA 

32 8 16 

a, + — a, + — , 

15/, ’ 3/, ^ 3/, 


•^10x11 ~-/^nxio 


EAC. 


EA 


K. 


EAC^ 


11x11 


EA 


23 


26 


11 


a, — a, — a, 

157 I 1. 


a, H a, H a, 

157 6 /, ' 3 / " 


(7.18) 
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Chapter 8 

Results and Discussion 

The first step in any aero-elastic response and stability analysis is the evaluation 
of natural frequencies of the rotor blade. Using the inertial and structural model 
developed in this study, a structural dynamic analysis was performed. It may be noted 
that the inertial and the structural operators given in Eqs.6.6 and 7 1 
respectively, are nonlinear. Since the structural dynamic analysis requires only linear 
terms, all the nonlinear terms are neglected The corresponding Enear equation for one 
beam finite element can be written as 


(8 1 ) 

Where [M]^ represents the mass matrix of /'* element, (given in Eq.6 6 and in 
Sec.6.2.1). The stififtiess matrix consists of three components. They arej^.^T'^J 
(given in Sec.6.2.3), (given in Eq. 7.5) and (given in Eq. 7.18). The 

element stiffness matrix is given as; 


The element matrices are assembled to form the global finite element model for 
the rotor blade. For the tip element, the corresponding local to global transformation 
(Ref. 8) is performed to take in account the sweep and anhedral angles. The 

transformation matrix [A"] is given in Appendix A. Imposing the root boundary 

conditions, the corresponding rows and columns from the global model are eliminated. 

The resulting matrix equation can written as; 
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(8.3) 


Performing an eigen analysis, the natural frequencies of an un-damped rotating blade in 
vacuum can be evaluated 

8.1 Validation 

In order to validate the finite element blade model developed in this study, the 
results of the present analysis are compared (for certain specific cases) with those 
available in the literature The data shown in Table 8.1 correspond to a uniform, 
untwisted straight blade for which uncoupled natural frequencies are available in the 
literature (Refs 1,3 and 6) The terminology soft-in-plate blade configuration indicates 
that the first non-dimensional rotating lag frequency is less than 1. 

Using the data given Table 8 1, the uncoupled natural frequencies of the rotating 
blade are calculated Natural frequencies are compared with those available in the 
literature (Refs l,3,and 6) as shown in Table 8.2. The results of the present analysis are 
in excellent agreement with those available in literature. The natural frequencies that are 
given in the Table 8 2 is for N=20. 


8.2 Practical data from the Industry 

The code is designed to work with non-dimensional quantities. Hence, the 
practical data is made non-dimensional. For that, mass per unit length (mo) is taken as 
8.45kg/m. Radius of the rotor blade (1) is taken as 6.6 m. Speed of the rotor(ao) ‘s taken 

as 314 rpm. 
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The non-dimensionalization has been done in the following manner 


El = GJ = 


El 


EA = 


EA 



— m 

m- — 
Wo 


The practical data obtained from the Industry contains non-uniform properties 
along the radial direction, and is given from 0 1136 (a non-dimensional value), as shown 
in Table 8.5. The variations in mass, mass moment of inertia and stif&iess are shown in 
the Figs. 8. 1-8.6. As a first attempt, two assumptions are made, i) the blade is fixed at 
0.1136, and ii) no variation of properties occurs within the elements. The results are 
tabulated in Table 8.6 In the next attempt, a linear variation of properties within the 
elements is considered. The results for this case are tabulated in Table 8 6. The results 
and mode shapes hence generated are as shown in figs. 8.7-8.10. 

The practical data also contains the root stiffiiess matrices from 0.0321 to 0 1136 
(non-dimensional), for flap, lag and torsional modes, implying that the blade is fixed at 
0.0321 (non-dimensional). The results hence simulated, are compared with the above- 
obtained results and are tabulated in Table 8.8. The effect of offset on frequencies is 
evident from the comparison. The mode shapes hence generated are as shown in figs. 
8.11-8.14 


Root stiffiiess matrix of the blade in lead-lag direction 


430.962 

-24.37 

-430 962 

-7.69 

-24.37 

2.042 

24.37 

-0.214 

-430.962 

24 37 

430^962 

7.69 

-7.69 

-0.214 

7.69 

0.796 
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Root stiffness matrix of the blade in flap direction 


141 78 
7 7102 
-141 78 
3 8683 


7.7102 
0 466 
-7.7102 
0.1629 


-141 78 
-7 7102 
141 78 
-3.8683 


3 8683 
0 1629 
-3 8683 
0 153 


Root torsional stiffiiess matrix 



6.951 
-8 056 
1.105 


-8 056 1.105 

9.756 -1 698 

-1.698 0.595 


8.3 Effect of root offset and Geometric pitch 

The effect of root offset (ei) is analyzed for the uniform blade. The frequencies are 
proportionally increasing with the root offset. A convergence study is made by giving a 
high stiffness from 0 to 0.0321 (non-dimensional), to simulate a root offset. The results 

converged at high stiffness [eI = 57680 j as shown in the Table 8 4. 

The variation of geometric pitch is taken as 12° at 0.25 and 4° at the end tip (i.e , 
non-dimension length 1). The coupled lead-lag and flap modes are observed. The lead-lag 
and flap frequencies reduced due to the geometric pitch. The coupled modes are shown in 
the Figs 8.15 and 8.16. 
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Table 8.1 Input data for uniform soft-in-plane rotor blade 


Uniform soft-in-plane rotor blade data 

- - - - ( non dimensional DroDerties'* 


0.0004 


0.000 


0 

m 

1 

Ps 

0 


0 

p. 

0 

0; 

0 

GJ 

0.001473 

EA 

20.0 

ei 

0.0 

^2 

0.0 

a. 

0.0 

Co = EACo/EA 

0.00021036 

EI,, 

0.0301 

EInn 

0.0106 


Table 8.2 Natural frequencies of uniform soft-in-plane rotor blade 


Mode 

Uncoupled natural frequency 

No. of elements 
(N = 20) 

Ref. [2] 

Ref.[l] 

Ref.[3] 

Ref.[6] 

l’‘Lag 

0.7311 

0.7311 

0.7326 

0.7317 

0.732 

2"“ Lag 

4.4530 

4.4531 

4.4563 

4.4825 


S'** Lag 

11.2868 

- 

- 

- 

- 

r‘Flap 

1.1244 

1.1251 

1.1247 

1.1245 

1.125 

2"'^ Flap 

3.4073 

3.4266 

3.4089 

3.4073 

- 

3"* Flap 

7.6171 

7.7154 

7.6376 

7.6218 

- 

1®* Torsion 

3.2633 

3.2633 

3.2632 

i 

3.263 

1®* Axial 

6.9389 

6.9389 

6.9533 

- 

- 












Radius 

3 21E-02 
1.14E-01 
1 21E-01 
1 .28E-01 
1.35E-01 
1.42E-01 
1.48E-01 
1.56E-01 
1.64E-01 
1.71E-01 
1.79E-01 
1.86E-01 
1.94E-01 
2.02E-01 
2.09E-01 
2.17E-01 


2.24E-01 


2.32E-01 


2.47E-01 


2.95E-01 


3.26E-01 


3.56E-01 


3.86E-01 



7.50E-01 


8.00E-01 


8.33E-01 

8.79E-01 

9.17E-01 

9.24E-01 

9.32E-01 

9.55E-01 

9.77E-01 

1.00E+00 


Table 8, 
M 

2 46E+00 
1.44E+00 
1.35E+00 
1 .27E+00 
1.21E+00 
1.14E+00 
1.07E+00 
1.02E+00 
9.93E-01 
9.61 E-01 
9.28E-01 
8.97E-C1 
8.66E-01 
8 85E-01 
9.05E-01 
9.24E-01 


9.43E-01 


5.98E-01 


6.29E-01 




6.89E-01 


1 .03E+00 


1.02E+00 


1 .02E+00 


1.02E+00 


1.02E+00 


1.02E+00 


1.02E+00 


1.01E+00 


1.01E+00 


1.01E+00 


1.01E+00 


1.00E+00 


1.00E+00 


1.00E+00 


9.99E-01 


9.96E-01 


9.95E-01 


9.94E-01 


9.64E-01 

9.27E-01 

8.95E-01 

8.54E-01 

8.14E-01 

5.94E-01 

4.18E-01 

2.49E-01 


.5 Industrial 

r 



2 12E-05 ~ 

5 73E-06 

5 41E-06 

5.19E-06 

5.76E-06 

6.66E-Q6 

7.58E-06 

8 59E-06 

8.94E-06 

8.86E-06 

8.78E-06 

8.67E-06 

8.56E-06 

8.45E-06 

8 29E-06 

8.10E-06 


7.91 E-06 


)ractical data of a helicopter blade 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


7.34E-06 


6.38E-06 

5.08E-06 

3.99E-06 

3.78E-06 

3.48E-06 

2.50E-06 

1.22E-06 

2.72E-07 


6 82E-04 
6.41 E-05 
5.32E-05 
4 40E-05 
3.91 E-05 
S.53E-05 
3.12E-05 
2 70E-0S 
2.50E-05 
2.44E-05 
2 38E-05 
2 31 E-05 
2.27E-05 
2.23E-05 
5.49E-05 
1.10E-04 


1.65E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.23E-04 


3.18E-04 
3.12E-04 
3.10E-04 
3.07E-04 
2 72E-04 
1.68E-04 
7.72E-05 
2.72E-06 


0 001990078 
0 001863175 
0 001892016 
0.002220812 
0.002549607 
0.002884171 
0.003189893 
0 003143747 
0.003091831 
0.003034148 


0.023189 

0.019093 

0.016094 

0.015805 

0.01194 

0.010095 

0.008306 

0.008018 

0.00773 

0.007441 


0.002979 

0.002863 

0.002703 

0.00265 

0.002689 

0.002825 

0.002975 

0.003047 

0.003059 

0.00305 


0.002976465 

0.007153 

0.003039 

0.002901476 

0.006864 

0.003029 

0.002838024 

0.008249 

0.003021 

0.002878403 

0.024227 

0.003107 

0.002918781 

0.040263 

0.00329 

0.00295916 

0.056241 

0.003591 

0.003016843 

0.079603 

0.003892 

0.003057221 

0.094601 

0.004479 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.003057221 


0.094601 


0.094601 


0.094601 


0.094601 


0.004672 


0.004672 


0.004672 


0.004672 




0.094601 


0.094601 


0.094601 


0.094601 


0.094601 


0.094601 


0.004672 


0.004672 


0.004672 


0.004672 


0.004672 


0.004672 


0.003057221 

0.094601 

0.004672 

0.003057221 

0.094601 

0.004672 

0.003057221 

0.094601 

0.004672 

0.003057221 

0.09^01 

0.004672 

0.003057221 

0.094601 

0.004672 

0.003057221 

0.094601 

0.004672 

0.003057221 

0.094601 

0.004672 

0.002861098 

0.094601 

0.004672 

0.002388094 

0.094024 

0.004255 


0.001626673 

0.001511306 

0.001222889 

0.001176742 

1.52E-03 

1.52E-03 


0.093447 

0.088832 

0.070374 

0.041994 

0.012517' 

0.012517 


0.002923 

0.002821 

0.002719 

0.001669 

0.000886 

0.002043 
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Table 8.6 Companson of natural frequencies for different assumptions in FE 
(ei= 0 0321, No. of elements = 30) 


Mode 

Linear variation within 
the elements 

No variation within the 
elements 

T' Lag 

1.0341 

1.0389 

2"“ Lag 

7.8315 

7.8479 

3^^* Lag 

21 4294 

21.461 

1^ Flap 

1.2296 

1.2043 

2"*^ Flap 

3.3896 

3.3176 

3^'* Flap 

6.3276 

6.2059 

Torsion 

5.6264 

5.5711 

1“^ Axial 

8.3336 

8.3967 


Table 8.7 Convergence of natural frequencies (ei = 0.0321) 


Mode 

No. of elements (N) =20 

No. of elements (N) =30 

l^Lag 

1.0343 

1.0341 

2"" Lag 

8.0399 

7.8315 

3^‘'Lag 

22 038 

21.4294 

1''‘ Flap 

12289 

1.2296 

2"^' Flap 

3.4134 

3.3896 

3"** Flap 

6.4059 

6.3276 

1*‘ Torsion 

5.6269 

5.6264 

Axial 

8.3421 

8.3336 


Table 8.8 Comparison of natural frequencies for different root offsets 


Mode 

Root offset ei = 0.0321 

Root offset ei = 0.1 136 

1** Lag 

1 0343 

1.269 

2"“ Lag 

8.0399 

8.295 

3"" Lag 

22.038 

22.823 

1** Flap 

1.2289 

1.435 

2"'* Flap 

3.4134 

4.092 

3*^^ Flap 

6.4059 

7.556 

1^ Torsion 

5.6269 

9.1577 

1®* Axial 

8.3421 

7.508 
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Figure 8.1 Distnbution of mass along the radius 
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Figure 8.2 Distribution of mass moment of inertia in flap direction along the radius 
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Figure 8.3 Distribution of mass moment of mertia in lead - lag direction along the radius 
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Figure 8.4 Distribution of flap wise bending stiffness along the radius 
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Figure 8.5 Distribution of lead - lag wise bending stiffness along the radius 
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Figure 8.6 Torsional stiffness distribution along the radius 
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Rgure 8.12 Mode shapes in flap(out of plane bending) mode 
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Figure 8.15 Coupled Lead lag mode 
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Figure 8.16 Coupled flap mode 
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Chapter 9 

Concluding Remarks 


A general structural dynamic model of a helicopter rotor has been formulated 
The model incorporates all the structural complexities present like pre-sweep angle, pre- 
twist, tip sweep angle, tip anhedral angle etc.. Due to practical requirements, in the 
present formulation, speed of rotation is taken as a variable. The rotor was modeled using 
beam finite elements having 14 degrees of fi'eedom The equations of motion are derived 
using Hamilton’s principle. The formulation was validated by comparing the results of 
the present analysis for a uniform hingeless rotor blade with that of those available in 
literature. Results have also been generated for practical helicopter rotor blade. 

The effect of root offset on the natural frequency and mode shapes of the 
hingeless rotor blade has been analyzed. The coupling effect of geometric pitch on the 
natural frequency and mode shapes of the practical rotor blade has also been analyzed. 

9.1 Scope for the Future Work 

Rotor blade aerodynamic model has to be formulated and integrated with the present 
structural dynamic model. The resulting analytical equations can be used for the study of 
aero-elastic response and stability of a rotor blade, coupled rotor-fuselage stability and 
vibration analysis of helicopters. 
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Appendix A 

Local to global coordinate transformation for tip element 


The local-to-global coordinate transformation for the swept-tip element can be writ- 
ten in the form 


{s°} = [A'•]{5^} {A.1) 

Where the subscript t denotes quantities associated with tip element; the superscripts 
L and G denote the local and global coordinate system, respectively; q is the vector 
of element nodal degrees of freedom, the transformation matrixs, [A^], is derived with 
constraint that the angular relationship between the swept-tip and the straight por- 
tion of the blade at the junction is preserved after deformation [Ref. 3 ]. For the 
translational degrees of freedom, the transformation is linear, as indicated by 



L 

f N 

1 “ 1 


1 ^ 

1 "" 1 

■ =%]' 


1 ^ J 

t 

1 ^ J 


(A.2) 


where the transformation matrix [Tg] is given by 
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(A.3) 


lTe] = 


cos Ag cos Aa — sin A, cos Ag sin A^ 
sin Ag cos Aa cos Ag sin Ag sin Aa 
- sin Aa 0 cos Aa 


An explicit form of the constraint relations for the rotational degrees of freedom 
is given by 


^ 1 

L 

r <!> 1 


> =irj' 

—to' 

V J 

t 

i I 


(A.4) 


Combining the Eqs. A.2-A.4 and rewriting in the form given in Eq. A.l, the ele- 
ments of transformation matrix are obtained. The elements of the transformation 
matrix are given below. 


Ill 



A^(l,l) 

= 

cos A, 

A^(l, 12) 

= 

sin A, cosAfl 

A^(2,9) 

= 

— sinAo 

A^(3,7) 

= 

sin A, sinAfl 

A^(4,4) 


cos Aa 

A^(5,5) 

= 

cosAa 

A^(6,2) 

= 

— sin A, sinAa 

A^(6,9) 

= 

— sinAj cosA^ 

A^(7,14) 

= 

— sin Ad 

A^(8,6) 

=z 

cos A, 

A^(9, 2) 

= 

cos A, sinAfl 

A^(9,9) 

= 

cos A, cosAfl 

A^(10,4) 

= 

5 cos A, sin Ad 

A^(10, 8) 

== 

1 sin Ad 

A^(ll,4) 

= 

cos Ad sin Ad 

A^(ll,ll) 

= 

cos Ad cos Ad 

A^(12,5) 

= 

cos Ad sin Ad 

A^(13,l) 

= 

- 5 sin Ad 

A^(13,5) 

= 

1 cos Ad sin A, 

A^(13, 13) 

= 

cos Ad cos Ad 

A^(14, 7) 

= 

cos Ad sin Ad 


A^(l,5) 

= 

sin Ad sin Ad 

A^(2,2) 

r= 

cos Ad 

A^(3,3) 

= 

cos Ad 

A^(3, 14) 


sin Ad cos Ad 

A^(4,ll) 

= 

— sin Aa 

A^(5, 12) 

= 

- sin Ad 

A^(6,6) 

= 

cos Ad 

A^(7,7) 

= 

cos Ad 

A^(8,4) 

= 

— sin Ad sin Ac 

A^(8, 11) 


— sin Ad cos A, 

A^(9,6) 

= 

sin Ad 

A^(10,2) 


1 cos Ad sin Ad 

A^(10,6) 

= 

1 sin Ad 

A^(10, 10) 

= 

cos Ad cos Ad 

A^(ll,8) 

== 

sin Ad 

A^(12,l) 


— sin Aj 

A^(12, 12) 

= 

cos Ad cos Ad 

A^(13,3) 

= 

-|sin As 

A^(13,7) 

zr 

1 cos Ad sin A, 

A^(14, 3) 

= 

— sin Ad 

A^(14, 14) 

= 

cos Ad cos Ad 
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Appendix B 

Modification of the Term [v,\] Associated with Kinetic Energy 


Expanding sm{6a + 4>k) and cos(0g + <f>k) , and assuming that ^k is small 

sin (0(5 + <}>k) ^ sin 6g + <^k cos Oq 
cos ( 0(3 + (pk) « cos 6a — <Pk sin 6 a 
Substituting the above approximationin [V^j], and rewriting as 

[^il ^ Wii] + [-^33 ]{^} 

where 

[■^ 33 ] = f cos^ 0(3 +lTn^ sin^0(j )(cos^0/ cos^Aa 

Jo 

— cos^ 0/ sin^ As sin^ Aa ) 

sin^ 0(5 + Irrij^ cos^ 0(5 )}{<t>q}{4>q}'^]dx 

= r [{(/mff — IfUrp, )sin0<5 COsOq COS^0/ COsAj COsAa 
(sin^ 0G — cos^ Oq ) cos^ 0/ cos As cos Aa 
+{Im^^ sin0G + Imrjr, cos0G)sine^ cosAo }{0g}lda: 
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